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1. Introduction 



< 

u 

^ ■ Relations between ellptic integrals and hyperelliptic integrals have been studied 

d . from the 19th century. For example Hermite [S] found the following formula 

(1.1) / , ''' =^[ '' . 

J ^y{z^-a){8z^-6az-b) 2^ J y/y^ ~ Say + b 

where y = {2z^ - b)/{3{z^ - a)). 

In this paper we derive several formulae which reduce hyperelliptic integrals to 
elliptic ones. They are obtained by comparing two expressions of global monodromies 
of the Heun equation. Here the Heun equation is the standard canonical form of a 
■^ . Fuchsian equation with four singularities, which is given by 

I: (1.2) ((f)\ (1 + ^ + ^) A + -?'"- ^ ) /» , 

j5 , \\dw J \w w — 1 w — tjdw w[w — l)[w — t)j 



with the condition 

(1.3) 7 + (5 + e = a + /3 + l. 



X 
^1 It is also known that solving Heun equation corresponds to the spectral problem for a 

certain model of quantum mechanics which is called the BCi Inozemtsev model. Set 

cuo = and uj2 = —Wi — u^. The BCi Inozemtsev model is a one-particle quantum 

mechanics model whose Hamiltonian is given as 

(1.4) F = -— + ^/,(/i + l)p(x + u;i), 

i=0 

where p{x) is the Weierstrass p-function with periods {2ui,2uJs), cuq, ^1,^2, 1^3 are 
half-periods, and k {i = 0,1,2,3) are coupling constants. This model is a one- 
particle version of the BCn Inozemtsev system jl] , which is known to be the universal 
quantum integrable system with S^v symmetry 0[7j. Let f{x) be the eigenvector of 
the Hamiltonian H whose eigenvalue is E, i.e. 

(1.5) (H - E)fix) = 1-^ + Y^ hih + l)p{x + iu,) - e) fix) = 0. 
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Then the couphng constants lo,li,l2j3 correpond to parameters a,/5, 7,5, e, a ratio 
of periods of elhptic function corresponds to a singular point t, an eigenvalue E 
corresponds to an accessory q, and an eigenfunction f{x) corresponds to a solution 
f{w). For details see 011111 [TI1[T2]. For the case Iq y^ and h = I2 = h = 0, Eq. (fT3|l 
is called the Lame equation. 

Hermite and Halphen investigated the Lame equation by hypothesizing that solu- 
tions are expressed by an elliptic Baker- Akhiezer (Block) type function and Krichever 
[S| described elliptic solutions to Kadomtset-Petviashvili equation by a Baker- Akhiezer 
type function. In this paper we investigate solutions to Eq. ()1.5p by considering the 
Hermite-Krichever Ansatz. In our situation, the Hermite-Krichever Ansatz asserts 
that the differential equation has solutions that are expressed as a finite series in 
the derivatives of an elliptic Baker- Akhiezer function, multiplied by an exponential 
function. More precisely, we are going to find solutions to Eq. (jl.5p of the form 



1^1.6) f{x) = exp {kx 



where $i(x, a) = exp{({a)x)a{x + tUj — a;)/cr(x + Ui) {i = 0, 1, 2, 3). 

Treibich and Verdier [1^ found and showed that, if k G Z>o for all i G {0, 1, 2, 3}, 
then the potential ^j^q h{h + ^)p{x+uji) satisfies the stationary higher KdV equation 
and they constructed a theory of elliptic soliton following Krichever's idea [Sj; while 
Gesztesy, Weikard [21 El, Smirnov P and the author [TUl[T2] obtained further results 
on this subject. Thus, the function ^j=o^«(^« + ^)pi^ + ^i) ^^ called the Treibich- 
Verdier potential; and is closely related with a hyperelliptic curve i/^ = —Q{E) where 
Q{E) is a polynomial in E which is determined for each Iq, li, I2, I3 G Z. For example, 
if the eigenvalue E satisfies Q{E) = 0, then the equation has a doubly periodic 
eigenfunction up to signs which corresponds to the Heun polynomial. (See [21 QQl 
ITT|.) In [12], global monodomies of Eq. (jl.5p with the condition Zq, ^1,^25 ^3 ^ ^ were 
calculated and they are expressed by hyperelliptic integral (see Proposition 13.11) . 

Belokolos, Eilbeck, Enolskii, Kostov and Smirnov studied the covering map from 
the hyperelliptic curve z/^ = —Q{E) to the ellitic curve p'{aY = 4:p{a)^ — g2p{a) — gs 
and they obtained relations among variables E,a and n in Eq. p.6|) for the case 
/o = 1,2,3,4,5, h = h = h = 0, the case /q = 2, /i = l,/2 = 0, /s = and the 
case Iq = 2,/i = l,/2 = 1,^3 = (see [Tj and the reference therein). By consider- 
ing the covering map they found transformation formulae like Eq. fjl.l|) that reduce 
hyperelliptic integrals to elliptic ones case by case. On the other hand, Maier found 
a pattern of the covering maps for the case of Lame equation (i.e. the case /q 7^ 0, 
^1 = ^2 = ^3 = 0) and conjectured formulae [UJ Conjecture L] by introducing the 
notions "twisted Lame polynomials" and "theta-twisted Lame polynomials" . 

In this paper we justify and develop the Hermite-Krichever Ansatz on the Heun 
equation without an advanced algebraic geometry technique for the case lo,h,h,h £ 
Z>o. Note that results on the Bethe Ansatz and monodromy formulae in terms of 
hyperelliptic integrals obtained in ^01 E] play important roles in our approach. As a 
result, the monodromies of Eq. (jl.5j) are expressed by elliptic integrals. To study the 
Heun equation by the Hermite-Krichever Ansatz, we need to consider the covering 
map in detail. For this purpose, we introduce twisted Heun and theta-twisted Heun 
polynomials that are based on Maier's ideas, and we obtain theorems that support 
the Maier's conjectures. By comparing two expressions of monodromies of Eq. ljl.Sp . 
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we establish transformation formulae between elliptic integrals of the first and second 
kinds, and hyperelliptic integrals of the first and second kinds. Hence, the mysteries 
of the elliptic-hyperelliptic integral formulae are unveiled by the monodromies. For 
the case /q = '^, h = h = h = 0, we obtain Eq. ljl.lj) as a transformation formula 
between elliptic integrals of first kind and hyperelliptic integrals of first kind. The 
formula for the second kind is written as 



, , I {2z'^ — a)dz 1 I'Sz^ — Qaz — b ^ f V^V 



^{z'^-a){Sz^-Qaz-h) SV z"^ - a 2VS J ^/y^ - 3ay + b' 

where y = {2z^ - b)/{3{z^ - a)). 

The Hermite-Krichever Ansatz would be applicable to the spectral problem of the 
BCi Inozemtsev model, because the monodromy is expressed in terms of an elliptic 
integral by applying the Hermite-Krichever Ansatz, and it is closely related with the 
boundary condition of the model. 

This paper is organized as follows. In section |2l we justify the Hermite-Krichever 
Ansatz on the Heun equation by applying results on the Bethe Ansatz and integral 
representation of solutions. In section 13 we obtain hyperelliptic-elliptic reduction 
formulae by comparing two expressions of monodromies. In section El we investi- 
gate zeros and poles of the covering map. In section we introduce twisted Heun 
and theta-twisted Heun polynomials, and obtain formulae which support Maier's 
conjectures. In section we give examples which cover the cases that the genus 
of the related hyperelliptic integral is less than or equal to three. Thus we obtain 
hyperelliptic-elliptic reduction formulae explicitly for more than 20 cases. 

Throughtout this paper we assume /q, h, h, h ^ ^>o and (/q, ^i, h, h) ¥" (0; 0; 0; 0)- 

2. Hermite-Krichever Ansatz 

In this section we review results on Bethe Ansatz and integral representation of 
solutions which are obtained in ^^, and apply them to justify the Hermite-Krichever 
Ansatz. 

Fix the eigenvalue E of the Hamiltonian H (see Eq. (ll.4j) ) and consider the second- 
order differential equation 

(2.1) {H - E)f{x) = (-^ + J2 k{k + l)p{x + uj,) - e\ fix) = 0. 

Let h(x) be the product of any pair of the solutions to Eq. ()2.1|) . Then the function 
h{x) satisfies the following third-order differential equation: 

(2.2) 

^-4 |^/,(/, + l)p(x + a;,)-^J ^-2 |^/,(/, + l)p'(x + a;,)J J h{x)=0. 

It is known that Eq. ()2.2|l has a nonzero doubly periodic solution for all E iili G Z>o 
(2 = 0,1,2,3). 

Proposition 2.1. TJl Proposition 3.5] If Iq,Ii,12,13 G Z>o, then equation i2.^} has 
a nonzero doubly periodic solution S(x, ii^), which has the expansion 

3 h-i 
(2.3) S(x, E) = co{E) +Y.Y. bf\E)pix + co^f-^ 

i=0 j=0 
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where the coefficients Co{E) and b^ (E) are polynomials in E, they do not have com- 
mon divisors and the polynomial cq{E) is monic. We set g = deg^ co{E) . Then the 
coefficients satisfy deg^ bj (E) < g for all i and j . 

We can derive the integral formula for the solution A(x, E) to Eq. (|2.ip in terms of 
the doubly periodic function E{x,E), which is obtained in [TU]. Set 

(2.4) 

Q{E) = E{x, Ef (^E-Y,^ hih + l)p{x + u;.)] + ls(x, E)^^^^ - \ {^^^ 

It is shown in jTUj that Q{E) is independent of x. Thus Q{E) is a monic polynomial 
in E of degree 2g + l, which follows from the expression for S(a;, E) given by Eq. ()2.3p . 
The following proposition on the integral representation of solutions is obtained in 

m-- 

Proposition 2.2. [10, Proposition 3.7] Let S(x,-E') be the doubly periodic function 
defined in Proposition \2. 1\ and Q{E) be the monic polynomial defined in Eq. /li2.4}) - 
Then the function 



(2.5) A(a:, E) = V^(^e^p J ^^.^^ 

is a solution to the differential equation h2.1\) . 

We consider the case Q{E) = 0. It follows from Eq. ()2.5|l that, A(a;, E)"^ = E{x, E). 
By considering zeros and poles, we obtain 



(2.6) A(x,E)^ = S(x,E) 



2 _ -/^ T7\ - •i^i=i 



{p{x) - eiY^{p{x) - e2Y^{p{x) - e^Y^ 



where /?i G {k, —k — 1} (z = 0, 1,2,3), C, ti, . . . are constants. From Eq. ()A.3j) . the 
function 

_ Yl^j^^^^+^-^P^^I^ a{x - t,)a{x + t,) 

^'^'^^ ^^'^' ^^ ~ a{xY^ai{xY'a2{xY'a:i{xY' ' 

is a solution to Eq.(j2IH). We set a = ^k=iPk^k for the case Q{E) = 0. Then we 
have 

(2.8) 

' A(x, E) (a = (mod 2cJiZ © 2CU3Z)) 



^^^ '^'" '' t exp(-2r7fca + 2tUfcC(a))A(x,E) (a ^ (mod 2cuiZ © 2CU3Z)) 

for fc = 1, 2, 3. 

Now we consider the case Q{E) 7^ 0. Then it is known that the functions A(x, E) 
and A(— X, E) are linearly independent. 

Set I = lo + h + h + h- The function S(x, E) is an even doubly periodic function 
which may have poles of degree 2/j at a; = a;i (i = 0, 1, 2, 3). Therefore we have the 
following expression; 

(2 9) -r. E) = ^o°^(^)n-^i(p(^)-Pfe)) 

^■' "^ ' ^ (p(3;)_ei)Mp(^)-e2)'^(p(x)- 63)^3 

for some values ti, . . . ,t/. It is shown in ^01 that, if Q{E) 7^ 0, then the values tj, 
—tj, Ui (mod 2u;iZ © 2c<j3Z) (j = 1, . . . , /, -i = 0, 1, 2, 3) are mutually distinct. Set 
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p{x) and Zj = p{tj). From Eq.(j23D we have (^^^) = ^§f^- We fix the 
? of tj by taking 

dE{x,E)\ _ 2^^Q{E) 



Let C(a^) be the Weierstrass zeta function, a{x) be the Weierstrass sigma function 
and crj(x) {i = 1,2,3) be the Weierstrass co-sigma functions. (See appendix.) If we 
put 2^J—Q{E)/'H.{x, E) into partial fractions, it is seen that 

2V-Q(g) 

It follows that 
f2.ll) 



A(x,E) 



\ rr3 w ^ ^. ^^P 9 Z^ Oogaitj + x) - loga(t 



x) - 2xC{t,)) 



bP{E)UU^(^ + ^^ 



exp I -x^Citj 



a{xy'^ai{xy^a2{xy^a3{xy^Ylj=i^itjJ \ i=i 

For the case Q{E) ^ 0, we set 



I 



(2.12) A(x, E) = iiz^L^ il — _^ y (^u 

a{xyoa,{xy^a,{xy^a,{xy-^UU^{t^) V tr 



Then the formula 



,(0)^ 



(2.13) A(x,E) = VC(^)A(x,i?) 

is derived. 

We establish the validity of the Hermite-Krichever Ansatz by using values ti, . . . ,ti. 
More precisely we show that, if lo,li,l2,h ^ ^>o, then an eigenfunction of the Hamil- 
tionian H with every eigenvalue E is expressed as in the form of Theorem 12.31 We set 
a = - Y!j=i tj + ELi ^kUJk and k = C(Ej=i tj - ELi ^fc^fc) - Ei=i C(^i) + ELi ^kVk, 
where rjk = C(^fc) (^ = 1? 2, 3). Since the set {—tj}j=i^„j is the set of zeros of A{x, E) 
and the position of the zeros and the poles are doubly-periodic, ±q; is determined 
uniquely mod 2uJiL © 2uJz1' and k is determined uniquely. From Eqs. ()A.51 IA.7|) it 
follows that 

(2.14) K{x + 2uJk,E) 

= exp I 2r7fe I ^ tj - ^ Ik'UJk' \ - 2ujk i^ ({tj) - ^ Zfc'^fc' | | A(a;, E) 
\ \j=i k'=i J \j=i k'=i J J 

= exp{—2rika + 2u!kC{ci) + 2KUJk)A{x, E) 
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for k = l,2, 3. We set 

(2.15) $.(x, a) = ^^"^ ,^ "^ ~ ""^ exp(C(a)x), (z = 0, 1, 2, 3). 

a{x + Ui) 

Then we have 

(2.16) l — j <^^ix + 2uJk,a) = exp{-2rika + 2uJkC{a))l — ] $i(x,a) 

for z = 0, 1, 2, 3, A; = 1, 2, 3 and j G Z>o. 

The following theorem asserts that eigenfunctions of the Hamitonian H are ex- 
pressed in the form of the Hermite-Krichever Ansatz. 

Theorem 2.3. Assume /q, ^i, h, h G Z>o and set I = lo + li + l2 + h. Let {—tj}j=i^__j 
be the set of zeros ofA{x, E) appeared in Eq.^EUB)- If ^ = ~ J2j=i '^i + Sfc=i ^^'^^ ^ 
{mod 2c<JiZ © 2UJ3Z), then we have 

(2.17) A{x,E)=exp{Kx) (^2^2^^ (/)' '^^^^' ") 

/or some values bj (i = 0, . . . , 3, j = 0, . . . , /j — 1). T/ie values a and k are expressed 
as 

(2.18) ^(") = ^' ^'^") = ^V^-^^' '^ = tM^^-^(^' 

where Pi{E), . . . , -P6(-E) are polynomials in E. 

If a = — X]i=i '^j + X]a:=i ^fc'^fc — {i^od 2ijj{L © 2ijjz£), then we have 

Cfc" 



(2.19) A(x,E)=exp(«:x)(c + X^5^6f(^yp(x + ^,) + E^^ f^ 

for some values c, Ck [k = 1, 2, 3) and bj (i = 0, . . . , 3, j = 0, . . . , Zj — 2) . 

Proof. Assume a ^ 0. From Eq. ()2.14|) and Eq. ()2.16p . the functions exp(Ka;) (^) $j(a;, a) 
(i = 0, 1, 2, 3, j G Z>o) and the function A{x, E) have the same periodicities. By- 
subtracting the functions exp(/€x) (^) $i(^)'^) from the function k{x,E) to erase 
poles, we obtain an holomorphic function that has the same periods as ^^{x^a) 
and it must be zero. Hence we obtain the expression ()2.17|) . For the case a = 0, 
we change {t^Y^^^ by ti -^ ti + a and tj -^ tj (j = 2,...,/). Then we can set 

a = - Y!j^i tj + Yl =i h^k = 0. By setting R = - Y!j=i C{tj) + ELi ^fc% we obtain 
the expression ()2.19|) . 

Next we investigate the values p(a), p'{a) and k for the case a ^ 0. The functions 

p(Ej=i ij - ELi ^kUJk). p\Y!j=i tj - ELi ^kUJk) and (k =)C(Ej=i ^i - ELi ^^^fc) " 
Ej=i C(^i) + Efc=i ^kVk are doubly-periodic in variables ti, . . . , t;. Hence by applying 
addition formulae of elliptic functions and considering the parity of functions p(x), 
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p\x) and C(3;)) we obtain the expression 

/ 3 

(2.20) p(5^t, - 5^4^fc) = Y. fnLjJ^(^^)^ • • • ' P(^i))p'(^n) ■ ■ ■ p'ih) 



J;/' 



1=1 fc = i h<h<-<Jm 

m: even 



I 3 

^'E^^- - E^'^^fc) = Y. fnLjJp(t^)^ • • • ' P(^O)p'fe) • • • P'it 

j = l fc = l Jl<J2<---<im 

m: odd 

Z 3 Z 3 

j=l fc=l j=l fc=l 

= E /il.Jp(^i)'---'P(^o)p'(t,j---p'(t 



ji/' 



3l)i 



m: odd 



c(fc) 



where fj^'j^{xi, . . . ,xi) {k = 1,2,3) are rational functions in xi,...x;. By ap- 
plying Eq. ()2.10p . the function p'{tj)/ ^J—Q{E) is expressed as a rational function 
in E and p{tj). Hence p(Ej=i^i - Hl=ihuJk), p'{Y!j=itj - Yl=ih^k) / \/ -Q{,E) 
and {C{Y!j^^tj - Yll=ih^k) - H]=iCitj) + YX=ihVk) / \/ -Q{E) are expressed as 
rational functions in the variable p(ti), . . . , p(t/) and E. They are symmetric in 
p(ti),...,p(t/). 

From Eq. (|2.9p we have 

(2.21) hf{E) \[{p{x) - Pit,)) = E{x, E){p{x) - e,y^{p{x) - e,)\p{x) - 63)'^ 

Hence the elementary symmetric functions Ylih<---<j , Pi'^h) ■ ■ ■ Pi'^ji/) (^' = 1; • • • ;0 

are expressed as rational functions in E. Therefore p{J2j=i '^j~J2k=i ^k^k), p'{^j=i tj — 

ELi ^^.O/v^^QM and (C(E5=i tj~ELi huJk)-EU C(^.)+ELi kVk)lV^Q{E) 
are expressed as rational functions in E. D 

Proposition 2.4. We have a -^ {mod 2ijj{L © 2ijj{£) as E ^ 00. 

Proof. We define the rational function S(z,i?) to satisfy 

(2.22) E{p{x),E) = E{x,E)/E^. 
From Proposition 12. 11 we have 

(2.23) E{z, E) = l + I,- -P(^.V;g) 

^ ' ^ ' E(z-ei)'i(z-e2)'2(^- 63)^3 

where P{x,y) is a polynomial in x and y. Let k G {1,2,3}. If E is sufficiently large 

and \z - Cfcl = \E\~^^^ , then we have \E{z,E) - 1| < 1^1"^/^ and \-§^E{z,E)\ < 
\E\~^^'^. Hence we have 



(2.24) ^ / 






Eiz,E) 







as \E\ -^ 00. Since the l.h.s. of Eq. ()2.24|) is equal to the number of zeros inside 

1 1 

1^ — Cfcl < \E\ 2{ifc+2) niinus the number of poles inside \z — e^l < \E\ ^{ifc+a) ^ jf |_£;| jg 

sufficiently large, then the function S(z, E) has Ik zeros in variable z inside the circle 
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k~efc| = \E\ 2(ife+2)^ j^ follows from Proposition l2.1l that the degree of P(2;, 1/E) in z 

is lo + h + h + h- Hence it is shown similarly that, if \E\ is sufficiently large, then the 

1 
function S(z, E) has Iq zeros in variable z outside the circle \z\ = |£'|2('o+2) . Therefore, 

if \E\ -^ oo, then 4 zeros of the function S(z, E) in z tend to Ck {k = 1, 2, 3) and /q 

zeros tend to infinity. 

Since the set {p(— ij)}j=i,...,Zo+ii+i2+«3 coincide with the set of zeros of E{z,E) in 

z, li elements of {— tj}j=i,...,io+ii+;2+«3 tends to Ui {i = 0,1,2,3). Hence we have 



a 



^(o+(i+(2+<3(_^^.) ^ ELi^fc^fc -^ (mod 2uJiZ © 2CU3Z) as E ^ 00. Thus we 
obtain the proposition. D 

3. Hyperelliptic-elliptic reduction formulae 

We obtain hyperelliptic-elliptic reduction formulae by comparing two expressions 
of monodromies. One is Eq. ()2.14|) . and the other one is the monodromy formula in 
terms of hyperelliptic integral that was obtained in 



Proposition 3.1. [12^ Theorem 3.7] Assume U G Z>o (i = 0,1,2,3^. Let k G 
{1,2,3}, qk G {0,1} andEo be the eigenvalue such that A{x+2uji;,Eq) = {—1)'^'^A{x,Eq). 
Then 

(3.1) A(a; + 2uJk, E) = (-l)^^A(x, E) exp — / ■^°+" ^ ' — rfE 

\ 2 jso ^ /-g(^) 



with e a constant so determined as to avoid passing through the poles in the integra- 
tion. 

Since the function p{x + Ui)"^ is written as a linear combination of the functions 



(^) ^("^ "*" ^*) ii ~ ^1 ■ ■ ■ 1^)1 the function S(x, E) can be expressed 



as 



(3.2) S(x,E) = c(E) + X:E4^^(^)(^) V(^ + ^. 



Set 

(3.3) a{E) = J24'{E). 

i=0 

From Prop osit ion 13 . 1 1 we have 

(3.4) M- + 2u„E) = (-l)-A(x,E)exp f-i f' -^V^E) + 2.,ciE) ^ 



Set 



(3.5) ^ = -2a - I ^^ ' dE, 

Ieo J-Q(E) 

(3.6) B = 2Cia) + 2K+ f ^^^' dE. 

Ieo J-Q{e) 
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By comparing Eq. ()2.14|) and Eq. ()3.4|) for the cases k = 1 and /c = 3, we have 

(3.7) r]iA + uJiB = 7iV^{qi + 2ni), 

(3.8) r/3A + cu35 = 7rv^(g3 + 2n3), 

for some integers rii and n^. By the Legendre's relation rjiu^ — rj^Ui = 71^/^/2 (see 
Ea. (|X3|l 1. we have 

(3.9) A/2 = (gi + 2ni)cu3 - {qs + 2n^)uJi, 

(3.10) - B/2 = (gi + 2ni)r/3 - (ga + 2^3)//!. 
Hence 

(3.11) "+0/ , dE = -{qi + 2ni)cu3 + (93 + 2^3)071, 



-Q(i5;) 



1 /'■^ c(F\ 
(3.12) C(«) + ^+0 I dE = -(gi + 2ni)r73 + (53 + 2^3)7/1. 

^^So J-Q{E) 



We set ^ = p(q;). By Proposition 12 .41 and the relation J{l/p'{a))dC, = J da, we have 

/•« de 1 /"^ a(E) ~ 

(3.13) / , = " = -9 / ^==rf^ 



Note that Q{E) is a polynomial with degree 2g + 1, and a{E) is a polynomial with 
degree g. Hence Eq. fj3.13p represents a formula which reduces a hyperelliptic integral 
of the first kind to an elliptic integral of the first kind. The transformation of variables 
is given by { = Pi{E)/P2{E) for polynomials Pi{E) and P2{E) (see Ea. dTT^ l. In 
sections |^ and |S1 we investigate the transformation of variables in detail. Let ao 
be the value a evaluated aX E = Eq. It follows from Eqs. (|3.5l 13. 9p that a^ = 
-(gi + 2ni)uJs + {q^ + 2113)101 and 

(3.14) a-ao + - / ^ dE = 0. 



If ao = (mod 2luiZ © 2u;3Z), then ({a - ao) = C(") + (^i + 2^21)773 - (^3 + 2n3)r]i, 
which follows from ({a + 2uJk) = ({a) + 27]}. {k = 1,3). Combining with Eqs. (j3.12l 
I3.14|) we have 

1 /-^ c(E) ~ (l [^ a(E) , 

(3.15) ^ = -n r^^ ^^ + ^o / r-^ dE 

^ -^^0 ^-Q{E) V ■^''' \J-QiE) 

If ao ^ (mod 2uo{L © 2CJ3Z), then C(«o) = —{Qi + 2ni)r]-i + (53 + 2n-i)r]i and 

(3.16) '^ = -^ / -^^^dE - C(a) + C(ao) 



-Q{E) 



■/i^o , /_g(i^) ^K) . /4|3 _ ^^1 _ ^3 
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Note that Q{E) is a polynomial with degree 2g + 1, and c{E) is a polynomial with 
degree g + 1. Hence Eq. (j3.16|) represents a formula which reduces a hyperelliptic 
integral of the second kind to an elliptic integral of the second kind. The following 
proposition presents asymptotic behaviour of p{a) and k as E -^ oo. 

Proposition 3.2. AsE-^oo, we have a ~ ^^ ^.Lo Wi+^), p(«) ~ -4^/(Eto ^*(^i+ 
1))2 and K ~ V^E{1 - 2/(ELo Wi + I)))- 

Proof. Set E{x,E) = J2i=o fii^)^^- ^Y definition we have fg{x) = 1. Substituting 
into ()2.2|) and comparing coefficients of E^, we have f'i{x) = v'{x)/2, where v{x) = 
'^i=Qh{k + ^)p{x + uJi). Hence /g_i(a;) = f(a;)/2 + C for some constants C. Therefore 
we have c{E) ~ E^, a{E) ~ E9-^Y.l=ySi + l)/2 as E ^ oo. From Eq.lQ we 
have Q{E) ~ E'^^+i as E ^ oo. Thus 



(3.17) 






c{E) ~_ [^ E^ 



(3.18) / / ^ ciE ~ / ^ — rfE ~ -2V^ 

as E —i' oo. By combining Eq. ()3.17|) with Eq. ()3.13|) we have a ~ ^ }^ J2i=(} ^ii^i + 1) 
as E ^ oo. It is known that p{x) -^ x~^ and ({x) -^ x~^ as x — ^ 0. Hence we have 
p{a) ~ -AE/iJ^l^hil, + 1))2 and C(a) ~ 2^/^E / {J^l^ l,{k + 1)) as E ^ oo. By 
combining with Eqs.dSHlEIISl) and Ea. (prTH|l we have k ~ ^/^E{1 - 2/(^^0 H^i + 
1)))- □ 

4. Covering map 

Let T be the elliptic curve whose basis periods are (2uji,2uj3) and Vq(^e) be the 
hyperelliptic curve defined by i/^ = —Q{E). In Theorem 12.31 a map n : Vq(£;) — > T 
which is called covering or covering map in PJ |B] is defined. Let a G T be the image 
of {E, I/) e Vq(b), i.e. a = 7i{E, v). From Theorem 12 . 31 and Proposition 13.21 we obtain 
the following expression of the covering map: 

^Nk{E) 

where Nk{E) and Dk{E) are monic polynomials such that deg^; Nk{E) = l+deg^; Dk{E). 
In this section we investigate zeros of the polynomials Nk{E) and Dk{E) {k = 1, 2, 3). 
The polynomials Nk[E) and Dk[E) [k = 1, 2, 3) can be calculated as the proof of 
Theorem 12.31 but practically it is hard to do. In P, the value p{a) is calculated by 
expanding functions as 

1 g2x'^ 

(4.3) -^--^— = -- -p{a)x+-p'{a)x'^ + ..., 

a{—a) X 2 D 

substituting them into Eqs. ()2.HI2.17j) and comparing coefficients of x-' (j = —2, —1, . . . 
althought it is not still effective. Maier found a pattern of the covering maps for the 



(4.1) p(a)-efc = - ^_3 ^^^ '\^^^^ ^^^ , (A; =1,2,3), 



(4.2) p{x) = - + ^ + 
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case of Lame equation and conjectured formulae of the covering maps [HI Conjec- 
ture L] by introducing the notions "twisted Lame polynomials" and "theta-twisted 
Lame polynomials". In this paper we calculate the polynomials Nk{E) and Dk{E) 
(A; = 1, 2, 3) by developing Maier's ideas. 

Fix the value E. Then the function K{x,E) and its zeros {— tj}j=i / are deter- 
mined, although we have ambiguity of choosing K{—x,E) and its zeros {tj}j=i,...,i- 
Hence the value a{= — X]j=i ^j + Tl,k=i ^k^k) (mod 2uJiL © 2u^Ij) is determined up 
to the sign ±. Note that this ambiguity corresponds to choosing v or —v on the 
hyperelliptic curve v'^ = —Q{E). The value p{a) is determined uniquely. 

Let p be an element of {0, 1,2,3} such that J2i=o^i^i — ^p (mod 2ti;iZ © 2ti;3Z). 
Assume that a = cOp (mod 2ljiZ © 2ci;3Z). Let pk{x) {k = 1,2,3) be cop function 
defined in Eq. ()A.7|) . It follows from Theorem 12 . 31 and Eqs. ()A.5UA.8|) that the function 
A{x,E) exp{—Rx)pi{xy^p2{xy^p3{xy^ is doubly-periodic. Hence the eigenfunction 
A{x, E) is written as 
(4.4) 

(Ei=o «i (^ - 62)^') + ^{z-ei){z-e2){z-e^) \Y.% bj {z - ea)^' 
A(x, E) = exp (kx) — T 



{z - eiY^/^iz - e2y^/^{z - esY^/^ 

where z = p{x) and /*^°^ (resp. l^^^) is the greatest integer that is less than or equal 
to (/o + li + h + k)/2 (resp. (/q + k + l2 + h- 3)/2). Let Ap (resp. Bp) be the set 
of eigenvalues E of the operator H such that the eigenfunction A(x, E) is written as 
the form ()4.4|) and satisfies k = (resp. k 7^ 0). 

Next assume a = Upi (mod 2ujiL © 2uj^'L) for p' G {0, 1, 2, 3} \ {p}. Let A; be an 
element of {1, 2, 3} such that uOpi ^ uOp + ujk (mod 2uJiL © 2ct;3Z). It follows from The- 
orem l2.3l and Eqs. ()A.5UA.8|) that the function A(a;, i?) exp(— fi;a;)pi(a:;)'ip2(a^)'^p3(a^)'^ 
has the same periodicity as the functions pk{x) and pki{x)pkii{x), where k',k" are 
determined as {A;, k' , k"} = {1, 2, 3}. Hence the eigenfunction A(x, i?) is written as 
(4.5) 

Vz - Cfc (E'=o S(^ - ^kY) + \/{z-ek'){z-ek") (Ei=o ^i(^ " ^fc^ 
A(x, i?) = exp (/tx) (^_e^)^i/2(^_e2);,/2(^_e3y3/2 

where 2; = p{x) and / (resp. /) is the greatest integer that is less than or equal to 
(/o + ^1 + ^2 + ^3 - l)/2 (resp. (/q + h + l2 + h- 2)/2). Let Ap> (resp. Sp/) be the set 
of eigenvalues E of the operator H such that the eigenfunction A(x, -E) is written as 
the form ()4.5|) and satisfies k = (resp. k 7^ 0). 
For the sets Aj, 5j {i = 0, 1, 2, 3), we have 

Proposition 4.1. (i) The eight sets Ai, Bi {i = 0, 1,2,3) are pairwise distinct. 

(a) The set Uf^^Ai coincides with the set of zeros of the polynomial Q{E). 

(Hi) Assume that Iq + h + I2 + h is even (resp. odd) and E e Bp. Then we have 

a[(o) 7^ (resp. bfio) ^ 0). 

(iv) Assume that Iq + h + h + h is even (resp. odd) and E G (Uf^giJi) \ Bp. Then we 

have bj 7^ (resp. a^ 7^ 0/ 

Proof. It follows from Theorem 12.31 that the values k^ and p{a) are determined 

uniquely for each E. li E & Ai U Bi, then we have a = uji (mod 2u;iZ © 2u;3Z) 
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and p{a) = < / ■ _ i 9 o\ Thus the sets Ai, Bi {i = 0,1, 2, 3) are separated by 

[^ Cj [I — i, z, oj. 

the values of k^ and p{a). Hence (i) is proved. 

Let JF be the space spanned by meromorphic doubly periodic functions up to signs, 

namely 

(4.6) -F= J-,,,,3, 

(4.7) J^ei,e3 = {f{x): meromorphic |/(x + 2u;i) = ei/(x), f{x + 2U3) = esfi^x)}. 

Then it is shown in [TUl Theorem 3.8] that Eq. ()2.1|) has a non-zero solution in the 
space J-" if and only if the value E' satisfies the equation Q{E') = 0. Suppose E' G Ai 
for some i. Then the value k or k is equal to zero. It follows from Eq. ()4.4|) or Eq. ()4.5p . 
that there exists a non-zero eigenfunction A(a;, E') satisfying A(x, E') e ^ej.eg for some 
ei, es e {-1, 1}. Hence we have Q{E') = and U^^^Ai C {E \ Q{E) = 0}.' 

Conversely suppose that Q{E') = 0. Then we can take the eigenvalue Eq in Propo- 
sition E^ as Eq = E'. From Eq. (j3.1ip we have a G cuiZ © uj^Tj. Hence the function 
A{x,E') is expressed as Eq. fl4.4|) or Eq. ()4.5|) . If k 7^ or /? 7^ 0, then the functions 
A(a;, E') and A(— x, E') are a basis of solutions to Eq. ()2.H) and they are not doubly- 
periodic up to signs. It contradicts the existence of a solution in the space JF. Thus 
we have k = or k = 0. Therefore we have U^^qAj D {E \ Q{E) = 0} and obtain (ii). 

We show (iii). Assume that a^^o) = (resp. 6|-(o) = 0). Then the degree of the pole 
at X = of the functions A(x, E) and A{—x, E) are less than Zq- Since R ^ 0, the 
functions A(x, E) and A{—x, E) are linearly independent. Hence the degree of the 
pole at x = of all solutions to Eq. (j2.H) must be less than Iq. But it contradicts that 
the exponents of solutions to Eq. (j2.H) are — /q and Iq + 1. Hence we obtain (iii). (iv) 
is proved similarly. D 

It follows from definitions oi Ai, Bi (i = 0, 1, 2, 3) that, if i? G Uf^g^j U -^ii then a 
solution to Eq. ()2.1|) is written as Eq. ()4.4|) or Eq. ()4.5|) . Conversely, it is shown that, if 
a solution to Eq. (P?T|) is written as Eq. (jO|) or Eq. fjO)) . then we have E G Uf^oAjUi^i. 
More precisely. 

Proposition 4.2. (i) If a solution to Eq. \2. 1]) is written as Eq.^^j, then we have 
E e ApU Bp. 

(ii) Let k G {1,2,3} and p' be the element of {0,1,2,3} such that ujpi = Up + Uk 
(mod 2c(JiZ © 2uj3'Z). If a solution to Eq. \2. 1]) is written as Eq. \4.5\j , then we have 
E e ApiU Bp, . 

Proof. We show (ii). Let A."'{x,E) be a solution to Eq. ()2.1|) written as Eq. ()4.5|) . Let 
Se be the set of all solutions to Eq. (j2.1|) for a fixed eigenvalue E. Then Se is a two 
dimensional vector space. Let M^/ {k' = 1, 3) be the transformations obtained by 
the analytic continuation x ^ x + 2ujki. Then the actions of My preserve the space 
Se and it is shown in ^21 §3.2] that the transformations Mk' do not depend on the 
choice of paths and M1M3 = M3M1. It follows from Eqs. (jXTl El) that 



(4.8) Mk'i'k'ix, E)) = exp{-2rik,ujp, + 2uk''k)k'{x, E) 

for some value k and k' = 1,3. 

Assume that Q{E) 7^ 0. Then it is shown in ^U] the functions A(x, E) and A{—x, E) 
defined in Eq. ()2.12j) form a basis of the space Se- Since the function A{x, E) is written 
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as Eq. ()2.14p . we have 

(4.9) Mfc,(A(x, E)) = exp{-2rik'a + 2iUk'X)Hx, E) 

for some value x ^i-nci fc' = 1, 3. It follows from changing variable x — > —x — 2uk' that 

(4.10) MkiM-x, E)) = exp{2r]k'a - 2uJk'X)H~x, E). 

From the assumption Q{E) ^ and pil| Theorem 3.8] we have Se^F = {0}. Hence 
we obtain that Mi or M3 have eigenvalues of neither 1 nor —1. Let M^ be the one 
that has an eigenvalue of neither 1 nor —1. Then the eigenvalues of M^ are distinct. 
On the other hand the function A"(x, E) is also an eigenvector of the operators Mj, 
(see Eq. fj4.8|) ). Hence the function A"(x, E) is written as CA(x, E) or CA(— a;, E) for 
some non-zero constant C. From Eqs. ()4.8M4.9|) or Eqs. ()4.8M4.10j) we have 

(4.11) ± {-2r]ia + 2uJix) = -2riiujpi + 2ujik + 27rv/^r2i 

(4.12) ± {-2ri^a + 2uj^x) = -'^Vs^p' + 2CU3K + 2^r^/^n3, 
for integers rii, 77.3. It follows that 

(4.13) (±a - ujp,){-2r]iuj-i + 2r]-iUJi) = 27r^/^{niUJs - UsUJi). 

From the Legendre's relation rjiuj^ — rj^uji = 7r-\/— 1/2 (see Eq. ()A.3|) ) and the relation 
ujpi = ±ujp> (mod 2co'iZ © 2u;3Z) we have a = ujpi (mod 2uJi'L © 2ci;3Z). Therefore 

E e Bp,. 

Assume that Q{E) = 0. Then the set of the eigenvalues of the operators M^/ 
{k' = 1, 3) on the space Se is either {1} or { — 1}, and the function A(x, E) in Eq. (j2.7p 
is an eigenvector of the operators M^/ {k' = 1, 3) with an eigenvalue given in Eq. ()2.8p . 
Hence the eigenvalue of the operator M^/ admits the expression as in Eq. ()4.9|) . On 
the other hand the function A"'{x,E) is also an eigenvector of the operators Mk' 
{k' = 1, 3) with an eigenvalue given in Eq. ()4.8|) . By comparing eigenvalues we obtain 
Eas. ffOnH:!^ and that a = Up' (mod 2cJiZ © 2CU3Z). Therefore E e Ap,. 

(i) is proved similarly. D 

The value p{a) — Ck {k = 1, 2, 3) is described by using the sets Aj, Bi {i = 0,k). 

Proposition 4.3. Let m{E') he the integer such that a{E) = (E - E')'^^^">a{E) 
and d{E') ^ 0, and m{E') be the integer such that Q{E) = {E - E')^^'^'^Q{E) and 
Q{E') y^O . We set 

(4.14) D{E)= JJ (^ _ ^')2+2m(s')--(i^') J] (E - E')'^""^'''^ 

E'eAo E'eBo 

(4.15) NkiE) = n (^ ~ E'f+2m{E')-ME') Y[{E- E')2+2m(i?')^ (fc = 1^ 2, 3). 

E'(^Ak E'eBk 

Then we have 

(4.16) p{a) = e, - 3 ^^^^^^ 

and deg^ Nk{E) = 1 + deg^ D{E) for k = 1,2,3. 

Proof. It follows from Theorem 12.31 and Proposition 13.21 that p{a) — Ck {k = 1,2,3) 
is expressed as Eq. ()4.1|l such that Dk{E) and Nk{E) are nionic polynomials, coprime 
and deg^ Nk{E) = 1 + deg^ Dk{E). We show that polynomials Nk{E) and Dk{E) 
are expressed as Eqs. ()4. 14114. ISj) . 
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We consider poles and zeros of the rational function p{a)—ek in variable E. Assume 
that a ^ 0,uJk (mod 2ujiZ(B2uj3'Z) . Then p{a) — ek is non-zero and non-infinity. Hence 
the eigenvalue E' such that a ^0,ujk does not satisfy D{E') = nor Nk{E') = 0. 

Next assume that E' is a zero of the rational function p{a) — Ck in variable E. We 
write p{a)-ek = {E-E')'^R{E) (m e Z>i, R{E') ^ 0, oo). Then E' i s an e igenvalue 
of the Hamiltonian that satisfies a = cjfc (mod 2uJiL®'2ujz'L) . From Eq. ()3.13p it follows 
that 

(4.17) / ^ = = — / / ^ dE. 

As E -^ E' , we have 

(4.18) r' '^ r. VPi-)-^^ = ciE - Er'\ 
Jek \/4:^^ - 92^ - 93 v(efe -efc/)(efe-efc//) 

(4.19) -I- r -^^3=dE = r^E - ET^'''^~'^^'^'^'''R{E)dE 
^Je' J-Q(E) J^' 



r. C2(E - E'V^""^^'^-"^^^'^/^ 



for some functions R{E) and non-zero constants Ci,C2- Hence we have m = 2 + 
2m{E') - m{E'). It follows from Proposition EH] that, if E' e -B^, then m{E') = 0. 
Therefore, if £" G Ak (resp. E' G B^), then the multiplicity of zeros of the polynomial 
Nk{E) at E = E' is2 + 2m{E') - m{E') (resp. 2 + 2m{E')). 

We finally assume that E' is a pole of the rational function p{a) — Ck in variable 
E. Then we have a ^ (mod 2ujiL © 2uj^'L) as E ^ E' . We write p{a) — Ck = 
R{E)/{E - E')"" (m e Z>i, R{E') ^ 0, oo). From Eq.(jSII3D it follows that 

(4.20) /"'"' 'g = -1 /' ,°'^) -^ 



-Q{E) 



As E —^ E' we have 



(4.21) / ^ =~ ^_ = Ci(E-EO'"/^ 
Joo V4^^ - 5-2^ - 6-3 V p(a) 

(4.22) - - /"" -^ELrfE ~ C'2(E - i^')i+-(i^')--(i?')/2^ 
2 ■^^' • / -Q{E) 



for non-zero constants (7i,C'2. Hence we have m = 2 + 2m{E') — m{E'). It follows 
from Proposition 14.11 that, if E' G Bq, then m{E') = 0. Therefore, if E' G Aq 
(resp. -E' G -Bq), then the multiplicity of zeros of the polynomial D{E) at E = E' is 
2 + 2m{E') - m{E') (resp. 2 + 2m{E')). U 

As a corollary we have 
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Corollary 4.4. Assume that a{E') ^ for all E' E U?^o(^* ^ ^i) '^'^'^ ^^^ equation 
Q{E) = does not have multiple roots. We set 

HiE- E') HiE- E')\ 

B'gAo E'&Bo 

- \[{E-E') n {E-E')\ (fc = 1,2,3). 

4Nk{E) 



(4.23) 




D{E): 


(4.24) 




Nk{E) 


Then we 


have 




(4.25) 







p{a) = Ck- 



(Elohik + lWDiE) 
and deg^ Nk{E) = 1 + deg^ D{E) for k = 1,2,3. 

Now we consider the function k. By combining Eq. ()2.18|) and Proposition 13.21 we 
obtain the expression 

2 \PJE)_ 



(4-26) /^= I 1- ^3 ,., . .. \-^^.V^Q{E) 



where Pk.{E) and Pk{E) are monic polynomials such that deg^; Pn{E) = deg^ P^{E) — 
g. We investigate zeros of the polynomial Pi^{E). 

Proposition 4.5. Assume that a{E') ^ for all E' G [Jf^Q{AiU Bi) and the equation 
Q{E) = does not have multiple roots. Then we have 

for some monic polynomials Pk{E) such that deg^ Pk{E) = deg^{Y[E'eA iE—E')Y[E'eB i^- 

E'))-g. 

Proof. It follows from Eqs. ()3.15l I3.16p . Proposition 14.11 and Corollary 14.41 that, if 
E' ^ Uf^Q(Aj U Bi), then the function k is holomorphic in E around E = E'. We 
investigate the degree of the poles at E' G Uf^Q^Ai U Bi). 

We show that, if fc G {1, 2, 3} and E' G Bk, then the function k is holomorphic in 
E around E = E'. From Corollarv 14.41 we have p{a) -^ Ck for some k G {1,2,3}, 
a ^ Uk and that ({a) is holomorphic around E = E'. It follows from Eq. ()3.12|l that 
the function k, is holomorphic in E around E = E'. 

Assume E' G Bq. From Corollarv 14.41 we have p{a) — Ck = gk{E)/{E — E'^ [k = 
1,2,3) for some rational function gk{E) such that gk{E') ^ 0, cxo. Since p{x) ~ 1/x^ 
and ({x) ~ 1/x as x ~ 0, we have a ~ (i? — E')/C' and C(tt) ~ C' /{E — E') for some 
non-zero constants C as E ^>- E' . From Eq. ()3.12|) . the function k has a single pole 
in E around E = E' . 

Next we consider the case E' G Aq. From Corollarv 14. 41 we have p{a) ~ {C'Y / {E — 
E') for some non-zero constants C as E -^ E' . Then we have C('^) ~ C!' /{E — E'Y^'^ 
as E ^ E' . On the other hand we have 

(4.28) / -A^)^dE^ [ R{E){E - E'Y^l^dE r^ C{E - E'Y'^ 

for some functions R{E) and constants C. Therefore we have k ~ —C'/{E — E'Y^'^ 
fromEa. dTT^ . 
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For the case E' E A^ {k E {1,2,3}), we have p{a) -^ e^ and C,{a) is bounded 
around E = E\ and f^, , ' dE is bounded around E = E' . Hence k is also 

•^^0 ^/~Q{E) 

bounded around E ^ E'. 

By combining with the fact that k, is expressed as Eq. ()4.2f)j) and that the equa- 
tion Q{E) = does not have multiple roots, the denominator of Eq. ()4.27|l must be 
expressed as Ue'^Ao {E - E') IIb'gBo (^ ~ ^') ' ° 

5. Twisted Heun and theta-twisted Heun polynomials 

To investigate Eqs. ()4.1(i^ I4.27J) and sets Ai, Bi {i = 0, 1,2,3), we introduce Heun, 
twisted Heun and theta-twisted Heun polynomials. 

To study the set A^ {i = 0,1,2,3), we consider Heun polynomials 0, which are 
related with quasi-exact solvability. For details see [TTJ §5]. For integers Po, Pi, P2, 
Ps such that —J2i=of^i ^ 2Z>o, Vgo,/3i,/32,/33 is defined as the vector space spanned 
by functions {pi(x)^ip2(a;)^'p3(a;)*p(a;)"}^^Q_ _^3_ ^^^3 ^^* "* ^ {-kJi + 1} (i = 
0,1,2,3) and 

•'1— aoil— Qijl— "2,1— 03) 2-/i=0 '^*/ '^>2) 

{0}, otherwise. 

If I0 + I1 + I2 + h is even, then the Hamiltonian H (see Eq. p.4|) ) preserves the spaces 

and if Iq + li + I2 + I3 is odd, then the Hamiltonian H preserves the spaces 

(5.3) f/_;„ _i^ _i2,/3 + l' ^-«0 -il,i2 + l-«3' f^-«0,'l+l,-«2 ,-«3' ^^^0 + 1 -«1 ,-«2 -«3 • 

Proposition 5.1. Let i G {0, 1,2,3}. Then there exists only one space ^ao,ai,o2,o3 ^'^ 
Eg. 115. ^} or Eg. 115. ,'3J) sttc/i t/ia^ X]fc=i c^fc^fc — ^j (^JTiod 2uji'Z(B2uj3'Ij). Moveover the set 
Ai coincides with the set of eigenvalues of the Hamiltonian on the space Uao,ai,a2,a-i- 

Proof. Assume that /o + ^i + ^2 + ^3 is even (resp. odd). Let f/ao,ai,a2,«3 ^^^ Ua' ,a' ,a', ,a' 
be different spaces in Eq. ()5.2|l (resp. Eq. ()5.3|l ). Then it follows directly that J2k=i ^k^k ^ 
Ylik=i ^'k^k (mod 2ti;iZ © 2u;3Z). Hence we obtain the first statement. It is shown in 
the proof of 1(7, Theorem 3.2] that there are no common eigenvalues on the spaces 

'^ ao ,cii ,a2 ,013 anQ L/ Q,^ Q,'^ Q,^ Q,^ . 



It is shown in jTHj that, if Iq + h + I2 + h is even (resp. odd), then the direct sum 
of the spaces in Eq. ()5.2p (resp. Eq. ()5.3p ). which we denote by V, is the maximum 
finite-dimensional subspace of JF, and the set of eigenvalues of the Hamiltonian on the 
space V coincides with the set of zeros of the polynomial Q{E). By combining with 
Proposition 14. H it follows that the set U^^qAj' coincides with the set of eigenvalues 
of the Hamiltonian on the space V. Let E & Ai. Then there exists a non-zero eigen- 
function A(x, E) written as Eq. ()4.4|l or Eq. ()4.5j) which satisfies R, = or k = 0. Then 

the function A(x, E) has the same periodicity as the function < \ , ) . ^ k o\ 

^ ^ \ Pi{x), (« = 1,2,3). 

Let UaoMi,a2,a3 be the space in Eq. ()5.2|l (resp. Eq. ()5.3|l ) such that J2k=i^k^k = ^i 
(mod 2uJiZQ)2iJ2,Z). Then the function /(x) G f/ao,ai,a2,a3 has the same periodicity as 

the function < '. ^ ,:. . k o\ Thus an eigenfunction in V^ whose eigenvalue is 
1 Pi{x), [t = 1,2,3). 



HEUN EQUATION IV 17 

an element of Ai belongs to the space f/ao,ai,a2,a3- Therefore we obtain that the set Ai 
coincides with the set of eigenvalues of the Hamiltonian on the space Uao,ai,a2,as- D 

For i G {0,1,2,3}, we denote the monic characteristic polynomial of H on the 
space Uao,ai,a2,a3 in Proposition 15.11 bv H^'^\E). We call them Heun polynomials. It 
follows from Proposition 15.11 that the set A^ {i = 0,1,2,3) coincides with the set of 
zeros of the Heun polynomial H^'^\E). 

Next we investigate the sets Bi (z = 0, 1, 2, 3) and introduce twisted Heun polyno- 
mials. Here we transform the Hamiltonian H. We set z = p{x), (f){x) = exp{Kx){z — 
eiY^^I'^{z - e2)"'2/2(2 - es)"'^/^, and R^, = 0(x)"^ o H o 0(x). Then we have 

(5.4) 5. ^ - 4(. - eO(. - e.)(. -e^^^t =^^} 

i=l 
+ {-{k + ll + l2 + h) {h +l2 + h-l0-l)z 

-K^ + eiik + h)^ + e2{h + hf + e3(/i + hf] 

j = l 

First we consider the set Bp. It follows from definition of Bp and Propositions 
14.11 14.21 that the condition E E Bp is equivalent to that there exists a solution to 
Eq. ()2.1|) written as Eq. ()4.4j) with the condition k ^ 0. We consider the situation that 
the function written as Eq. ()4.4|l with the condition R ^ satisfies the differential 
equation Eq. (j2.H) . We set 

m m 

(5.5) $(x,E) = Y^ aj{z - e^y + \/ {z - e-^){z - e2){z - e^)^hj{z - e^Y , 

j=0 j=Q 

ttj = (j > /(°\ j < 0) and bj (j > l^°\ j < 0), where 1^°^ (resp. 1^°^) is the greatest 
integer that is less than or equal to {Iq + I1 + I2 + h)/'^ (resp. (/q + ^1 + ^2 + ^3 — 3)/2). 
Then Eq.(|2IH) for /(x) = A{x,E) in Eq.(jl3I) is equivalent to 

(5.6) (z - e^){z - e2){z - 63) (^« - ^) ^(x, E) = 0. 

By substituting Eq. (j5.5j) into Eq. (j5.6|) and comparing coefficients of {z — €2)^^^^ and 
A/(^^-"e^J(^^^"e^y(F^^"e^(z — 62)^'^'^, we obtain relations 

(5.7) {h + k + h-lo-l- 2j)(/i + ^2 + ^3 + ^0 - 2j)aj 

i+4 j+3 

+ ^ aj'CaaiJ,f,K)+R ^ bj>CabiJ,f,K)=Q 

(5.8) ih + k + h-lo-^- 2j)ih + l2 + l3 + lo-^- 2j)b, 

j+4 j+4 

j'=j+l j'=j+2 

where CaaiJJ', k.), c^iUJ', z^), ci^iJJ', ^) and CuU,]', R) are polynomials in ei, 62, 63, 
R and E. It follows from Proposition 14. II that, if $(a;, E) satisfies Eq. (j5.6|) . R y^ and 
lo + h + h + h is even (resp. odd), then the value ai^o) (resp. bjio)) is non-zero. We 
set a|-(o) = 1 (resp. b^o) = 1) for normalization. 
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If Iq + h + h + h is even (resp. odd) and Iq > li + I2 + I3 — 4: (resp. /q > 
^1+^2+^3-1), then values 6;-(o)_2,af(o)_i,fef(o)_3, • • • , ao (resp. %o)+i, 6[(o)_i, %o), • • • , ao) 
are determined recursively by Eqs. ()5.7l I5.8|) . To satisfy Eq. ()5.6|) the coefficients of 
{z — 62)-' and \/{z — ei){z — e2){z — e3){z — 62)-' (j = 0, 1, 2, 3) on the l.h.s. should 
be zero. Hence we have simultaneous equations for k,E, which are Eqs. ()5.7|l5.8j) for 
the cases j = —1, —2, —3, —4 with a_4 = 6_4 = ■ ■ ■ = a_i = 6_i = 0. We solve them 
by throwing away the solutions corresponding to the case h = and removing the k 
factors by the elimination method in the theory of Grobner bases. Then we obtain 
an equation Ht^^\E) = 0, where Ht^^\E) is a polynomial which is normalized to 
be monic. We call Ht^\E) the twisted Heun polynomial. (See [0] for twisted Lame 
polynomial.) Note that we can remove completely the term including k, because it 
follows from Proposition 14.31 that the number of eigenvalues E whose eigenfunction 
is written as Eq. ()4.4|) is finite. 

Let us consider the case ^0+^1+^2+^3 is even (resp. odd) and Iq < /1+/2+/3— 4 (resp. 
^0 < ^1 + ^2 + ^3- !)• We set A = 6(ii+«2+i3-i„-4)/2 (resp. A = a^i^+i^+i^^ig^i)/2). Then 
values (ij and bj except for 6(;^+/2+;3_;o_4)/2 (resp. a(^i-^+i^^i^^iQ_iy2) are determined 
recursively by Eqs. (j5.71 15. 8p . To satisfy Eq. (|5.6|) the coefficients of {z — 62)^ and 

^/{z^^^'e^i){z^^^'e2)iz^^~'^{^ ~ ^2)^ {j = 0, 1, 2, 3) on the l.h.s. should be zero. Hence 
we have simultaneous equations for k, E^ A. We solve them by throwing away the 
solutions corresponding to the case k = and removing the R and A factors by the 
elimination method in the thoery of Grobner bases. Then we obtain an equation 
Ht^P\E) = 0, where Ht^^\E) is the twisted Heun polynomial which is normalized to 
be monic. Note that we can also completely remove the term including R and A. 

Note that elements of the sets Ap sometimes appear as solutions for the recursive 
equation for the case k = 0, although some elements in Ap might not appear as 
solutions because the condition a;(o) 7^ (or 6[{o) 7^ 0) may be broken. 

We consider the sets Bp' {p' G {0, 1,2,3} \ {p})- Let k be an element of {1,2,3} 
such that ujpi = cjp + cj^ (mod 2ti;iZ © 2uj3'L). It follows from the definition of Bpi and 
Propositions 14.11 14.21 that the condition E G Bpi is equivalent to that there exists a 
solution to Eq. (j2.ip written as Eq. (j4.5j) with the condition k 7^ 0. We consider the 
situation that the function written as Eq. (j4.5p with the condition k 7^ satisfies the 
differential equation Eq. ()2.H) . We set 
(5.9) 



^{x,E) = y/z-ek I ^aj{z-eky + y/{z - ek'){z - Ck") ^bj{ 
. j=o / \ i=o 



z - ekV 



^j = (j > ^5 j < 0) and bj (j > I, j < 0) where k', k" are determined as 
{fc, fc', k"} = {1,2, 3}, z = p{x) and / (resp. /) is the greatest integer that is less than 
or equal to (/q + h + h + h - l)/2 (resp. (/q + li+l2 + h- 2)/2). Then Eq.(|2ll} for 
f{x) = A{x,E) in Eq. ()4.5j) is equivalent to 



(5.10) {z - e,){z - e2){z - 63) (h. - E^ <l>(x, E) = 0, 
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(see Eq. fl5.4|) ). By substituting Eq. ()5.5|) into Eq. ()5.6|) and comparing coefficients of 

^ z — tk^z — CkY'^'^ and \/Jz'^-e^7}iz'^^'^^i^ ~ 6^)-'"'"^, we obtain relations 

(5.11) {k + k + h-lo-^- 2j)(/i + /2 + ^3 + ^0 - 1 - 2j)s 

j+4 j+4 

+ X^ a-fCaaiJ, f, k) + K ^ bj,CabiJ, f , K,) = 

(5.12) (/, + /2 + /s - /o - 3 - 2j){h + l2 + h + lo-2- 2j)bj 

j+4 i+3 

+ X] h'^bbij,j',fi) + K ^ ajvCba(j,j', k) = 0, 

i'=i+i i'=i+i 

where CaaiJJ', /«), CabiJJ', k), CbaiJJ', k) and Cbb{j,j', k) are polynomial s in c i, eg, eg, 
K and £". It follows from Proposition 14.11 that, if ^{x,E) satisfies Eq. ()5.10p . k 7^ 
and Iq + h + h + h is even (resp. odd), then the value bf (resp. aj-) is non-zero. We 
set bf = 1 (resp. a^ = 1) for normalization. 

If Iq + h + I2 + h is even (resp. odd) and Iq > li + I2 + h — 2 (resp. /q > 
^1 + ^2 + ^3 — 3), then values af_i, &[_i, af_2; • • • , cto (resp. 6^-, a.f_i, &[_i, . . . , Oq) are 
determined recursively by Eqs. ()5.111 I5.12|) . To satisfy Eq. ()5.10|) the coefficients of 
y/z — ek{z — CkY and a/(-2 — ek'){z — ek"){z — CkY (j = 0, 1, 2, 3) on the l.h.s. should 
be zero. Hence we have simultaneous equations for k, E. We solve them by throwing 
away the solutions corresponding to the case k = and removing the k factors by 
the elimination method in the theory of Grobner bases. Then we obtain an equation 
Ht^P\E) = 0, where Ht^^^E) is the twisted Heun polynomial which is normalized 
to be monic. Note that we can also completely remove the term including k. 

If IQ + I1 + I2 + h is even (resp. odd) and Iq < I1 + I2 + I3 — 2 (resp. /q < ^1 + ^2 + ^3 — 3) , 
then we need to set A = a{i^+i^+i^^i„_2)/2 (resp. A = b(j^+i^+i^_i^_3)/2) and obtain 
simultaneous equation for k, E^ A. The twisted Heun polynomial Ht^^ KE), which is 
normalized to be monic, is defined similarly. 

It follows from the definitions and Propositions 14. 11 14.21 that the set of zeros of the 
Heun polynomial H^'^'> (E) coincides with the set Ai and that the set of zeros of the 
twisted Heun polynomial Ht^^\E) coincide with the set Bi (i = 0, 1, 2, 3). 

The value p{a) is expressed by Heun and twisted Heun polynomials. 

Theorem 5.2. Assume that the Heun polynomials H^'^\E) and the twisted Heun 
polynomials Ht^'^\E) do not have multiple zeros and that their zeros do not intersect 
with zeros of a{E) for generic periods {2uji,2uj3). Then we have 

(5.13) ,(«) ^ e. - 3 ^^"^(^)^f"(^)^ 

ij:lohih + l)yH(oKE)Ht(o){Ey 

and <legEH^^\E)Ht^^\Ef = 1 + deg ^ H^°\E)Ht^°\Ef for k = 1,2,3. 

Proof. Assume that {2uji,2uj3) are the basic periods that satisfy the assumption of 
the theorem. Since the set Ai (resp. Bi) {i = 0, 1, 2, 3) coincides with the set of zeros 
of the polynomial H^^\E) (resp. Ht^^\E)), we obtain 

(5.14) H^'\E) = Yl {E - E'), Ht^'\E) = Y[ {E - E'). 
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From Corollary 14.41 we obtain Eq. ()5.13p for the case that satisfies the assumption of 
the theorem for basic periods. On the other hand p{a) is a rational function in E. 
Hence Eq. fl5.13j) is true without generic conditions. D 

In section ini it is shown that, ii g < 3, then the assumption of Theorem 15.21 is true. 
Now we introduce the theta- twisted Heun polynomial H6{E) to calculate the poly- 
nomial Pk,{E) appeared in Eq. ()4.27j) . 

We consider the case ^j=o ^i^i — ^ (j^od luj^L © luj-^T). We set 

$o(a;,a)(p(a;) - es)^ 



(5.15) u 

(5.16) V, 



' (p(x) - ei)'i/2(p(x) - e2)'^/2(p(x) - e3)'3/2 ' 
(^$o(a;,a))(p(x)-e2)^' 



(p(x) - ei)'i/2(p(x) - e2)'^/2(p(x) - 63)^3/2 ' 

(see Eq. (l2.15|l ). Let A(x,-E') be the function written as Eq. ()2.17|) and assume that 
K = and a ^ (mod cuiZScusZ). It follows from periodicities and position of poles 
that, if ^j=Q liOJi = (mod 2ijj{L © 2c<j3Z), k = and a ^ (mod ujxL © co'3Z), then 
the function A(x, E) is written as 

i I 

(5.17) A(x, E) = Y, c,u, + Y, djVj 

j=0 j=0 

where / (resp. /) is the greatest integer that is less than or equal to (/0+/1+/2+/3 — 1)/2 
(resp. (/o + ^1 + ^2 + ^3 — 2)/2). Conversely it is shown similarly to Proposition 14.21 
that, if a ^ (mod tUiZ © Ci;3Z) and the r.h.s. of Eq. ()5.17|) satisfies the differential 
equation ()2.1|) . then we have k = 0. 

We substitute Eq. (l5.17|l into Eq.dSHJ. By using a relation (^) $o(x, a) = {2p{x) + 
p{a))^o{x, a) and comparing coefficients of Mj+4 and Vj+4, we obtain relations 

(5.18) {h + l2 + h-lo-2- 2j){h + l2 + h + lo-l- 2j)cj 

j+4 j+A 

+ X] 9'Ccc(i,i',a) + p'(tt) Y dj'Ccd{j,f,a) = 

j'=j+l j'=j+2 

(5.19) ih + k + h-lo-S- 2j)(/i + l2 + h + lo-2- 2j)dj 

j+4 j+4 

+ Y dj>Cdd{j,f,a) + p'{a) Y^ Cj>Cdc{j,j',a) =0, 
j'=j+i j'=j+i 

where Ccc{j,f,a), Ccd{j,f,a), Cdc{j,f,a) and Cdd{j,j',a) are polynomials in d, 62, 
63, p{a),p'{a) and E. If A(x, £J) satisfies Eq. ()2.H) . a ^ (mod cJiZ © ci;3Z) and 
lo + h + h + h is even (resp. odd), then it follows that Q{E) 7^ and the value df 
(resp. Cj) is non-zero, which is obtained similarly to Proposition 14. II We set df = 1 
(resp. Cf = 1) for normalization. 

If Iq + li + I2 + I3 is even (resp. odd) and Iq > li + I2 + I3 — 2 (resp. Iq > li + I2 + 
Is — 3), then values q, di_^, q_;^, . . . , Cq (resp. d^^^, q-,^, (ij-_2, . . . , Cq) are determined 
recursively by Eqs. (j5. 18115. 19p . To satisfy Eq. (j2.1|) the coefficients of Mj and Vj (j = 
0,1,2,3) on the l.h.s. should be zero. Hence we have simultaneous equations for 
p{a), p'{a),E. Note that the relation p'(a)^ = 4(p(a) — ei){p{a) — e2){p{a) — 63) 
is also satisfied. We solve them by throwing away the solutions corresponding to the 
case a = (mod u;iZ©ct;3Z) and removing the p{<y), p'{<y) factors by the elimination 
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method in the thoery of Grobner bases. Then we obtain an equation H6{E) = 0, 
where H6{E) is a polynomial normalized to be monic. We call H6{E) the theta- 
twisted Heun polynomial, (see jHj for theta-twisted Lame polynomial.) Note that 
we can remove completely the term including p{a),p'{a), because it follows from 
Eq. ()2.18|) that the number of eigenvalues E whose eigenfunction is written as a linear 
combination of Uj and Vj (j = 0, 1, . . . ) is finite. 

Let us consider the case ^0+^1+^2+^3 is even (resp. odd) and Iq < /1+/2+/3— 2 (resp. 
lo<h + h + h-^)- We set A = qi^+i^+i^_i^_2)/2 (resp. A = d(^i^+i2+i^_i^_sy2). Then 
values Cj and dj except for C(/i+i2+«3-«o-2)/2 (resp. d(^i^+i2+i^_ig_sy2) are determined 
recursively by Eqs. (j5.18tl5.19|) . To satisfy Eq. (|2.ip the coefficients of Uj and Vj (j = 
0,1,2,3) on the l.h.s. should be zero. Hence we have simultaneous equations for 
p{a),p'{a),E,A. We solve them by throwing away the solutions corresponding to 
the case a = (mod UiZ © ci^sZ) and removing the p{a),p'{a),A factors by the 
elimination method in the theory of Grobner bases. Then we obtain an equation 
H6{E) = 0, where H6{E) is the theta-twisted Heun polynomial which is normalized 
to be monic. Note that we can also remove completely the term including p{a), p'{a), 
A. 

Next we consider the case X]i=o ^i^i — ^p (^^od 2ti;iZ © 2co'3Z) and p 7^ 0. We set 



/^ ^„x ^ , X aix — a — ijJr^ ( ( ^, X 1 p'(tt + i^D) A A 

5.20 vl; X, a) = ^ ^^-^ exp C a + ^p + 7: ^^ . '' x 

a[x) W 2 p{a + ujp) - ep J J 



Then the function \E'p(x, a) has the same periodicities as the function ^q{x, a){p{x) 
ei)~^^l'^[p[x) — e2)~'^/^(p(x) — 63)"'^/^, which follows from Eq. ()A.4|) for the case z 
X + Up and z = —ujp. We set 



(5.21) uj 



^p(x,a)(p(x) -62)^' 



(5.22) t; 



{p{x) - ei)'i/2(p(a;) - e2y^/^{p{x) - esY^/'^' 
{£^p{x,a)){p{x)-e2y 
' {p{x) - ei)'i/2(p(x) - e2y-/'{p{x) - 63)^3/2 ' 



Let A{x,E) be the function written as Eq. ()2.17p and assume that k, = and a ^ 
(mod uJiZ © uj^Ij). It follows from periodicities and position of poles that, if k = 
and a ^ (mod tuiZ © uj^Z), then the function A(x, E) is written as 



I I 

(5.23) ~A{x,E) = J2cju, + J2dj 

j=0 j=0 



Vj. 



Conversely it is shown similarly to Proposition 14.21 that, if a ^ (mod uJxL © ct;3Z) 
and the r.h.s. of Eq. ()5.23|) satisfies the differential equation ()2.1|) . then we have k = 0. 



22 KOUICHI TAKEMURA 

We substitute Eq. ()5.23|) into Eq. ()2.1|) . By using relations of elliptic functions and 
comparing coefficients of Mj+4 and f j+4, we obtain relations 

(5.24) {p{a + Up) - epfih + k + h - k - 2 - 2j){h + /s + ^3 + ^o - 1 - 2j)c,- 

j+4 i+4 

+ X^ Cj>CcciJ,j',a) + p{a + ujp){p{a + ujp)-ep) ^ dj^CcdiJ,j',a) = 

(5.25) {p{a + Up) - ep)\h + /a + ^3 - ^o - 3 - 2j){h + h + h + k - 2 - 2j)rf,-4 

i+4 i+4 

+ X] ^j"^dd{jj', a) + p'{a + ujp){p{a + Up) - Cp) ^ Cj'CddJJ', a) = 0, 
i'=i+i i'=i 

where CcdJJ', a), c^dUJ', a), CddJJ', a) and CddiJJ', a) are^polynomials in d, eg, eg, 
p(q; + Up), p'{a + Up) and E. It is shown similarly that, if A{x, E) satisfies Eq.(j2IH), 
a; ^ (mod tUiZ © t^gZ) and /q + ^1 + ^2 + ^3 is even (resp. odd), then the value di 
(resp. ci) is non-zero. We set di = {p{a + Up) — e^)^' (resp. q = {p{a + Up) — e^)^') 
for normalization. 

Similarly to the case Yli=o ^i^i — (iiiod 2co'iZ©2ci;3Z), the coefficients Cj and dj are 
determined and we have simultaneous equations for p{a + ujp), p'{a + ujp), E and so on 
to satisfy Eq. ()2.1|l . Note that, if /o+/i+/2+^3 is even (resp. odd) and Iq < I1+I2+I3 — 2 
(resp. /o < ^1+^2 + ^3-3), then we need to set A = C(i^+i^+i^_i^_2)/2 (resp. A = 
(i(;j+i2+i3_;Q_3)/2) and obtain simultaneous equations for p{a + Up), p'{a + Up), E, A, 
else we obtain simultaneous equations for p{a + u!p), p'{a + u!p), E. By throwing away 
the solution a = (mod tuiZ © co'3Z) and removing the factors except for E by the 
elimination method in the theory of Grobner bases, we obtain the theta-twisted Heun 
polynomial H6{E), which is normalized to be monic and the equation H6{E) = 0. 

Theorem 5.3. Assume that for generic periods {2ui,2uj3) the polynomials a{E), 
Q{E), He{E), H^\E), Ht^\E), {i = 0, 1,2,3) do not have multiple zeros, and their 
zeros do not intersect each other. If deg^ H9{E) — deg ^ H^^\E)Ht^^\E) = —g, then 
we have 

Proof. Assume that {2uji,2uo3) are the basic periods that satisfy the assumption of 
the theorem. From Proposition 14.51 and Theorem 15. 2t we have 



If the eigenvalue E' such that a ^ (mod cuiZ © uj^T) satisfies H9{E') = 0, then 
we have k, = and Pk{E') = 0. By assumption, the equation H6{E) = does not 
have multiple roots. Hence the polynomial Pk{E) is divisible by H6{E). On the 
other hand we have deg^ He{E) = -g + deg ^ H^^\E)Ht^^\E) = deg ^ Pn{E) from 
the assumption. Proposition 14.51 and Eq. ()5.14j) . Hence we have Pk{E) = H6{E) and 
Eq. ()5.2fij) for the case that satisfies the assumption for basic periods. On the other 



hand n/ ^J—Q{E) is a rational function in E. Hence Eq. ()5.26p is true without generic 
conditions. D 

In section ini we show that, if (7 < 3, then the assumption of Theorem 15.21 is true. 
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6. Examples 

In this section we explicitly calculate covering maps and several functions that 
have appeared in this paper. In section 16.11 we review constructions of covering 
maps and related functions. In section 16.21 we observe relations among different 
couphng constants (/q, ^i, ^2, ^s)- In sections 16.3116. 61 we express covering maps and 
related functions for all the cases where the genus of the hyperelliptic curve is less 
than or equal to three, and some cases of genus equal to four. By substituting them 
into the ones in section 16.11 we obtain several explicit formulae. Especially , for 
each case the monodromies are expressed as Eq. ()6.Hl and the transformation formula 
between elliptic integrals of the first kind (resp. second kind) and the hyperelliptic 
integrals of the first kind (resp. second kind) are expressed as Eq. (|6.4|) (resp. Eq. (|6.5|) 
or Eq.(jniSI))- 

6.1. Review of covering maps and related functions. We review constructions 
of the covering maps and related functions; later we explicitly express their form for 
each case. 

The doubly periodic function S(x,-E) is defined in Proposition 12. II Based on the 
function H(a;, ii^), the polynomials Q{E), c{E) and a{E) are determined in Eqs. ()2.4l 
13. 2113. 3|) . Let A(x, E) be the solutions to Eq. (j2.1|) which is expressed as Eq. (j2.12|) and 
Eq be an eigenvalue such that Q{Eq) = 0. Then there exists qi,q2,(l3 £ {0, 1} such 
that A(x + 2uJk, Eq) = (— l)''*A(x, Eq) for A; = 1, 2, 3. The monodromies of solutions 
to Eq. ()2.1|) are expressed by hyperelliptic integrals as 

(6.1) A(x + 2u„ E) = (-l)-A(x, E) exp { -\ [" z]:nM^+2^^dE 

\ ^-^^^ \J-Q{E) 

for k = 1,2,3 (see Eq.Q). The Heun polynomials i7»(E) {i = 0,1,2,3), the 
twisted Heun polynomials H&\E) [i = 0, 1,2,3) and the theta-twisted Heun poly- 
nomial H6{E) are introduced in sectional Let g be the genus of the hyperelliptic curve 
u"^ = —Q(E). From explicit expressions of functions a(E), c{E), H9{E), H^'^\E) and 
Ht^^'{E) {i = 0,1,2,3) given in a list in sections 16. 3116. 5| assumptions of Theorem 
15.21 and Theorem 15.31 are shown to be correct for the case g < 3. It can be shown 
similarly that they are also correct for the cases Iq < 10 and li = I2 = h = 0. In the 
case where assumptions of Theorem 15.21 and Theorem 15.31 are correct, the image of 
the covering map ^(= p{a)) and k are calculated as in Eqs. ()5.13[l5.26p . i.e., 

4iJ"W(E)iJ't(fe)(E)^ 

(E -=0 hih + l)yH(^KE)Ht^^) (Ef 



(6.2) C = efc o ^ ' , , ^ ^ (A; = 1,2,3), 



2 \ H9{E) 



Then a formula which reduces a hyperelliptic integral of the first kind to an elliptic 
integral of the first kind, i.e., 

(6.4) / , = -- / ±^ dE{= a) 

J 00 ^4|3 _ g^l _ ^3 ^Joo ^-Qi^E) 

is obtained (see Eq. ()3.13|) ). Next we express a formula which reduces a hyperelliptic 
integral of the second kind to an elliptic integral of the second kind. Since the value Eo 
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satisfies Q{Eq) = we have Eq G Ai for some i G {0, 1, 2, 3} (see Proposition l4.1|) . Let 
,^0 be tlie value ^ evaluated aX, E = Eq determined by Eq. (j6.2p . If Eq E AiU A2U A-^, 
then we also obtain a formula which reduces a hyperelliptic integral of the second 
kind to an elliptic integral of the second kind, i.e., 

(,5) lf<El=,E.-.,f^i=. 

(see Eq. (l3.1(ij) ). If Eq G Aq, then we have 

(6.6) '^ = -9/ r^^ dE + C\- / ±J— dE 

^ -^^0 ^-Q{E) V ■^''' \/-QiE) 

(see Eq. ()3.15|) ). Thus the solution to Eq. ()2.H) is expressed in the form of the Hermite- 
Krichever Ansatz (see Eq. fj2.17p ) and the monodromies are expressed as 

(6.7) A(x + 2iUk, E) = exp(-2r/fca + 2tUfcC(tt) + 2KuJk)A{x, E) 

for A; = 1,2,3 (see Eq. (|2.14j) ). where the value a is determined from Eq. (|6.4p . The 
value K is written as Eq. ()6.3|) . By expressing a and C(q^) ^-s elliptic integrals in variable 
^(= p(a)), the monodromies are represented in terms of elliptic integrals through a 
transformation by the covering map ^ (see Eq. fl6.2|) ). 



6.2. Relations among the cases of different coupling constants. We comment 
on relationship among several coupling constants {lo-ih-ihih)- Let Sj^'^^^'^^'^;^ (x, £") 

(^resp. ^ig^i.^i^^i^ [X), 0-10,1^,12,13 W, %,i^,i2,i3 {X}, ^,1^,1^,12,1^ , f^io,i^,i2,i3 ) oe tne lunctions 
E{x,E) (resp. Q{x), a{x), c(x), ^, k) for the periods (2a;i,2u;3) and the coupling 
constants {lo,h,h,h)- 

The following proposition is obtained by parallel transformation x ^ x + Uk {k = 
1, 2, 3) and definitions of functions (see Eqs.Q EH EH ElEH 1^^ ^: 

Proposition 6.1. We have 

(6.8) 

2;:i£f (X, E) = Hi^-Sf (X + u„E), Hg-Sf (X, E) = ^£S{x + u,, E), 

^!^:^{x.E)=^!>^:^{x+u,^Ei 

-j^(2(Ji,2w3)/ N _ -j^(2(^i,2u)3)/ x _ -jr^{2a;i,2aj3)/ N _ -y(2tjji,2u)z) , •. 

(2lji,2w3) _ {2a;i,2aJ3) _ (2tJi,2a;3) _ (2lji,2w3) 
^lo,h,h,h ~ ^h,lo,l3,h ~ ^h,l3,lQ,h ~ ^l3,h,h,lo ' 

for X = Q, a, c and x = (,, i^- 

By Propotision 16. II and changing periods, it follows that 
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Proposition 6.2. We have 

(6-9) 2!--j---)(x, E) = ^^-'-^'{x + u„ E) = ^l^^ix, E), 

2aS;^(x, E) = ^£:^ix + c,, E) = ^li:^ix, El 

T^(2aJi,-2<^i-2w3)/ \ _ Tr^(2<^i,-2a;i-2a;3)/ \ _ Tr^(2a;i,2w3)/ n 

'0)'3)'2,'l ^ ' '2,'1,'0.'3 ^ ' '0,'l,'2,'3 ^ '' 

-j^(-2<^i-2aJ3,2<^3)/ N _ -jr^(-2i^i-2<^3,2c<;3)/ N _ -jr^(2a;i,2aj3)/ n 
^lo,l2,li,h y-^) ~ ^h,h,l2,lo y-^) ~ ^lo,h,l2,h v"^^' 

{2uJi-2uJi-2uJ3) _ {2uJi-2uJi-2uJ:i) _ (2aJi,2aJ3) 
^lo,h,h,h ~ ^h,lo,l2,h ~ ^lo,h,h,h ' 

(2a;3,2tJi) ^ (2LJ3,2a;i) ^ (2a;i,2a;3) 
•^/0,^3,'2,il -^^2,^1,^0,^3 -^'0,'l,'2,'3 ' 

(-2a)i-2aJ3,2aJ3) _ (-2(^i-2aJ3,2aJ3) _ (2a;i,2a;3) 
^lo,l2,h,h ~ ^l3,hM,lo ~ ^lo,h,h,l3 ' 

for X = Q, a, c and x = ^, k. 

We can replace the coupling constants {lo,h,l2,h) to the case h > h > h > h 
by repeatedly applying Proposition lfi.2[ We consider examples mainly for the case 
^0 > ^1 > ^2 > ^3 in sections IU31I6. 61 

The relationship between the case (/q, 0,0,0) and the case {Iq,Iq,Iq,Iq) is con- 
nected by changing the scale of the periods (2ti;i, 2ci;3) ^^ (^1,^3), and the relation- 
ship between the cases (/o,/i,0,0) and (Zq, ^o; ^i; ^1) is connected by changing periods 
(2u;i,2a;3) ^ (u;3,2a;i). 

6.3. The case g = 1. If the genus of the related hyperelliptic curve z/^ = —Q{E) is 
one and the condition Iq > h > h ^ h is satisfied, then we have three cases 

(/o,/i,/2,/3) = (1,0,0,0), (1,1,0,0), (1,1,1,1). 

6.3.1. r/ie case (/o,^i,/2,/3) = (1,0,0,0). We have 

E{x,E)=E + p{x), a{E) = l, c{E) = E, Q{E) = {E + ei){E + e2){E + e^), 

i7(0)(E) = 1, if(^)(E) =E + ek, {k = 1, 2, 3), 

Ht'-'\E) = 1, (i = 0, 1, 2, 3), H9{E) = 1. 

For this case, an eigenfunction of the Hamiltonian H with eigenvalue E is expressed 

as 

exp(C(a;)x)cr(x — a)/a{x), E = —p{a), 

and we have k = 1 and ^ = —E. The formulae ()6.4[ 16. 5|) are trivial, because they 
only change signs by ^ = — -E. 

6.3.2. r/ie case (Zo,^i,^2,/3) = (1,1,0,0). We have 

E{x,E) =E + p{x) + p{x + ui)-3eu a{E) = 2, c{E) = E - Se^ 

i7(o)(E) = E-4ei, H'^^^E) = E^ - 2eiE + g2 - llel H^^^E) = H'^^^E) = 1, 

Q{E) = H^'^\E)H^^\E), Ht^'^\E) = m^^\E) = 1, 

i7t(2) ^E)=E + 562 + 363, Ht^^'> {E) = E + 362 + 563, He{E) = 1. 

The formulae ()6.4[I6.6|1 are related with Landen transformation, i.e. changing periods 

(2u;i,2t^3) 4-^ (u;i,2t^3). 
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6.3.3. r/ie case (Zq, ^1, ^2, ^3) = (I5 15 15 !)■ We have 

E{X, E)= E + p{x) + p{x + Ui) + p{x + U2) + p(x + CJg), 

a(E)=4, c{E) = E, Q{E) = {E + 4ei){E + 4e2){E + 4es), 

H^^){E) = Q{E), i7W(E) = l, (A; = l,2,3), i/t^^^) ^ ^^ 

m(*=)(E) = ^2 + ScfcE + 4(72 - 32e2, (fc = 1, 2, 3), /f^(E) = E^ - ^g^. 

The formulae ()fj.4[l(i.(ij) are related with the changing scales of the periods {2uJi, 2uj^) *-^ 

(^1,^3)- 

6.4. The case g = 2. If the genus of the related hyperelliptic curve u'^ = —Q{E) is 
two and the condition /o > ^1 > ^2 > ^3 is satisfied, then we have seven cases 

r (2,0,0,0), (2,1,0,0), (2,1,1,0), (2,2,0,0), 
lio,ii,i2,^3J-| (2,2,1,1), (2,2,2,2), (1,1,1,0). 

6.4.1. The case [hMMM) = (2,0,0,0). 

E{x,E)=9p{xy + 3Ep{x) + E''-lg2, a{E) = 3E, c(E) = E^ - §^2, 
H(^\E) = E^-3g2, H^''\E) = E - 3ek, (A; =1,2,3), QiE) = UUh'^'KE) 
Ht^°\E) = l, Ht'^''\E) = E + 6ekAk = 1,2,3), H9{E) = 1. 

By setting ^ = —y/6, E = z, g2 = a/3, gs = 6/54, Eq. ()l.l|l is obtained from 
Eq. (j6.4|) . and Eq. (jl.7j) is obtained from Eq. (j6.5j) . 

6.4.2. The case {Iq, /i, /a, ^3) = (2, 1, 0, 0) . 

E{x, E) = 9p(x)2 + 3{E - 2ei)p{x) + {E + 4ei)p(x + Ui) + E^ - hciE - 21e\, 
a{E) =AE- 2ei, c{E) = E^ - he,E - 27e? + f ^2, Q{E) = [iLo H^'^iE), 
H^^\E) = E + 4ei, H^^\E) = E'^ + 2{Ae^ + 3e2)^ - 526263 - 126^ - 316^, 
H^^\E) = 1, H'^^^E) = E^ + 2(462 + 363)^ - 526263 - 126^ - 3l6i, 
m(°)(E) =E- f 61, Ht^^\E) = E^- eiE - 1016? + f (72, 
i/t(2) Ie)=E + 1462 + 563, Ht^^^ {E) = E + 562 + 1463, H9{E) = 1. 

6.4.3. The case (/q, ^1, ^2, ^3) = (2,1,1,0). 

E{x, E) = dpixf + 3{E + 2e3)p{x) + {E + ci - 8e2)p{x + ooi) 

+ {E + 62 - 86i)p(x + UJ2) + E^ + besE + 486^ - f ^2, 

a(E) = 5E + 1363, c(E) = ^2 ^ 563^ + 486^ - 15(72, QiE) = ULoH^'KE), 

H^^^{E) = 1, H^^^{E) = E -862 + 61, H^'^\E) = E + 62 - 861, 

H('\E) = E^ + 36sE^ + (516^ - 20^2)£^ + f 63(72 - ^^3, 

Hti^){E)=E^ + f6sE+'-f6l-'-fg2, 

Ht^^\E) =E^ + 2(-62 + 36i)E + 336^ + 746i62 - 1036?, 



Ht^'^\E) =E^ + 2(-6i + 362)^ + 336? + 746162 - 1036: 
Ht^^\E) =E + 1763, H9{E) = 1. 

6.4.4. The case {k, h, k, h) = (2, 2, 0, 0) . 
E{x, E) = 9p(x)2 + 3{E - 66i)p{x) + 9p{x + cji)^ + 3{E - 66i)p{x + uji) 

+E^ - IbeiE - 366? + 1^2, 
a{E) = 6{E - 661), c{E) = E^ - 15eiE - 366? + 6c/2, 

H^°\E) = E{E^ - 12eiE - 1446? + 12^2), H^'^\E) = E"^ - 18eiE - 276? + 9^2, 
H^'^^E) = 1, {k = 2, 3), QiE) = H'^^\E)HW(E), 
Ht^^^ {E)=E- I861, m(i) {E) = E^- 4326? + 36^2, 

Ht^^\E) =E^ + 9(562 + 363)^2 + 21661(61 - 62)^ + 432(61 - 62)(62 - 63)(63 - 61), 
ift(3)(^) =E^ + 9(362 + 563)^2 + 21661(61 - 63)^ - 432(61 - 62)(62 - 63)(63 - 61), 

H9{E) = E'- f 61E - f 6? + f (72. 
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6.4.5. The case [Iq, h, k, h) = (2, 2,1,1). 

E{x, E) = 9{p{xf + p{x + uJif) + 3{E - 4ei){p{x) + p{x + ui)) 

+ {E - lQei){p{x + UJ2) + p{x + 003)) + E^- lOeiE - Sle? - |^2, 
a{E) = 8{E-7ei), c{E) = E^ - WciE - 51el - 3g2, H^^\E) = E - 16eu 
if(i)(E) = E^-AeiE^ -2{7g2 + 57el)E^ -4:{105el + 31g3)E + 81gl-llUg3ei-127e\, 
Q{E) = H(°\E)H^^\E), H^^\E) = 1, H'-^^E) = 1, 
Ht^^\E) = {E + llei)(E2 - UciE + 49e? - 27^2), 
Ht^'^^E) =E^ + 22eiE - Sdlej + 81^2, 

Ht^^\E) = E^ + (1263 + 4462)^^ + (3726263 + 1746i - 186^)^2 
+28(62 + 63)(-19l62 - IO6263 + 736i)E 

+ (63132636i - 607726162 + 70626^61 - 274956^ + 344576^) 

Ht^^\E) = E^ + (1262 + 4463)^^ + (3726263 + 1746^ - 1861)^2 
+28(62 + 63)(-19l62 - IO6263 + 736^)^ 

+ (631326261 - 607726^63 + 70626^61 - 274956^ + 344576^) 
He{E) = E^- fe,E - ^6? - fg2. 

6.4.6. The case {k, hM, h) = (2, 2, 2, 2) . 

S(x, E) = 9{p{xf + p{x + uj,f + p{x + uj2f + p(x + uj^f) 

+3E{p{x) + p{x + uji) + p{x + CU2) + p{x + UJ3)) + E^- 27^2, 
a{E) = 12E, c{E) = E"^ - 24^2, 

Q{E) = {E~ 126i)(E - 1262)(^ - 1263)(^2 _ 48^2), H^'^^E) = Q{E), 
Ht'~^\E) = {E + 2Aei){E + 24:62) iE + 24:es), H^^^E) = I, (A; = 1,2,3), 
ift«(E) = ^6 + 72ekE'' + 324(-862 + g2)E^ - 3456(^726^ + ^73)^' 

+31104(8^26^ - ^g^ek - gl)E' + 746496 (8 ^73^26^ + 4^71 + gl - Sgiel), 
He{E) = E^- '^g2E' + '^g,E^ + '-^glE^ - ^^^2^3^? - '-^{gl + 4^1)- 



6.4.7. The case (loJuhJs) = (1,1,1,0). 

E{x, E) = {E- 363)p(x) + {E- 3e2)p{x + c^i) + {E - 3ei)p{x + CJ2) + E^ - |^2, 

aiE) = 3E, c{E) = E^-\g2, Q(E) = [ito^^^H^)' 

H^'>\E) = E^-3g2, H<^^\E) = E-3ei, H^^\E) = E-3e2, H^^\E) = E-3e3, 

Ht^^\E) = l, Ht^^^E) = E + 6ei, Ht^^\E) = E + 6e2, Ht^^\E) = E + Ge^, 

Hd{E) = 1. 

Hence the functions except for H(x, E) are the same as the ones in the hst for the case 

(/q, ^1, h, h) = (2, 0, 0, 0). Especially the monodromies of the solutions to Eq. (j2.1|) for 

the case {lo,h,l2,h) = (2,0,0,0) and the one for the case (loJiyhJs) = (1,1,1,0) 

are exactly the same. For the case (Zq, h, h, h) = (0, 1, 1, 1), we have 

E{x, E) = {E- 3ei)p{x + uoi) + {E - 3e2)p{x + UJ2) + {E - 3e3)p{x + LU3) + E^ - §^2 

and the rest of the functions are the same as the ones in the list. 

6.5. The case g = 3. If the genus of the related hyperelliptic curve u"^ = —Q{E) is 
three and the condition /q > ^1 > ^2 > ^3 is satisfied, then we have 14 cases 

r (3,0,0,0), (3,1,0,0), (3,1,1,0), (3,1,1,1), (3,2,0,0), 
[lo.hMM = { (3,2,1,0), (3,2,2,1), (3,3,0,0), (3,3,1,1), (3,3,2,2), 
[ (3,3,3,3), (2,1,1,1), (2,2,1,0), (2,2,2,0). 

We omit some twisted Heun polynomials because we do not need to use them to 
express the covering map. 
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6.5.1. The case {I0J1J2J3) = (3,0,0,0). 

E{x, E) = 22hp{xf + AhEp{xf + Q{E^ - f ^2)^(0:) + E^ - Ibg^E - ^g^, 

a{E) = Q{E^ - f ^2), c{E) = E' - fg^E - if ^3, 

i7(o)(E) = E, H^^\E) = ^2 + QetE + 45e^ - 15^2, {k = 1, 2, 3), 

Q(S) = UUh^HE), Ht^'KE) =E^- f ^2, 

i7t(^)(E) = E^ + Ibe^E + f ((72 - 12e2), (A; = 1, 2, 3), H9{E) = 1. 

6.5.2. The case (/q, /i, /2, ^3) = (3,1, 0, 0) . 

E{x, E) = 225p{xf + 45{E - 2ei)p{xf + {6E^ - 24Eei - 201el - fg2)p{x) 

+ {E - 862 - 3e3){E - 362 - 863)^(2; + c^i) 

+^3 _ 7g^^2 _ io46?E - fg2E - 726? - ^g^. 
a{E) = 7E^ - UciE - 1526? - 10^2, 
c{E) =E^- 761^2 - (1046? + lg2)E - 1026? - if (73, 

QiE) = UUh^'KE), H(°KE) = E'- 16eiE - 326? + 5^72, 

E^\E) = ^3 - 961^2 - (1176? + 4^2)^ + 696? - 188^3, 

H^'^\E) = E - 362 - 863, H^'^\E) = E - 862 - 363, 

fft(o)(^) = ^2 _ ^g^^ _ 737g2 _ 100^^^ Ht^^\E) = E^ + lle^E - 6266? + 50(^2, 

H9{E) = E-17ei. 

6.5.3. The case (/q, h, h, h) = (3, 1, 1, 0). 
E{x, E) = 225p{xf + 45(263 + E)p{xf + §(8^2 + 32^63 + 3326^ - 175^2)pC 

+ {E^ + (861 - 1062)^ + 206? - 846162 - 796^)p(a; + tui) 
+ (E^ + f-106i + 86,)^ - 796? - 846169 + 20e'^)i3(x + UJ2) 



x) 



+ {E^ + (861 - 1U62JJ^ + 2U6f - 846162 - 796^jp(a; + 
+ {E^ + (-IO61 + 862)^ - 796? - 846162 + 206^)p(x 

+^3 + 7^2g^ + (i45g2 _ 59^^)^ ^ |g^^^ _ 5|5^^ 

a{E) = 8E^ + 26^63 + 1406^ - 85g2, Q{E) = Ilto^^'H^), 

c{E) =E^ + 763^2 + i45g2^ _ 2^g^E + A2esg2 - 255g3, 

H^^\E) = E^ + 1263^2 + 4(goei - 17g2)E + ?>2{2e^g2 - hg^), 

H^^\E) = E^ + (862 - 106i)E + 206^ - 846261 - 796?, 

H^^\E) = E^ + (861 - 1062)^ + 206? - 846162 - 796^, H'^^^E) = 1, 

Ht('\E) =E^ + f 63i? + ^6^ - ^^2, 

m(3) (E) = E^ + 4763^3 + (2436^ - f (72)^' + 80(7(^3 - 10(7263)^ 

+ 193326^^2 - 3125^1 - 33908^363, 
He{E) = E^ + fe,E + if^6i - f^(72. 

6.5.4. The case (/q, ^1, h, h) = (3, 1, 1, 1). 
E{x, E) = 225p{x)' + 45Ep{xf + 6{E' - 'ig2)p{x) + E^ - '-fg2E - '-^g, 

+ {E - 15e2){E - 15es)p{x + cui) + {E - l5ei){E - 15e^)p{x + 0^2) 

+ {E - l5ei){E - 1562)p(x + u^), 
a{E) = 9(^2 - 15(72), c{E) = E^ - 45(72^ - 540(73, Q{E) = UU H^'HE), 
/7(o)(^) = E^- 54c/2^^ - 864^3^ - 135^|, Ht^^\E) = E^ - 75^2, 
i7W(E)=E-156fc, (A; = 1,2,3), 
Ht^^\E) = E^ + ASckE^ + 54(^2 - 86^)^2 _ 8540e3E - 2025(3^1 ^ 48^gg^ _ 54^ 



/J^(E) = ^3 _ 53^2^ _ 540^, 



3- 



6.5.5. The case {lo,li,l2,h) = (3,2,0,0). 

S(x, E) = 225p{xf + 45(E - 66i)p(a;)2 + 3(2^^ - 24E6i - 1536?)p(a;) 
+9{E + 96i)p(a; + u^f + 3{E - 126i)(E + 96i)p(a; + cui) 
+E^ - 2161^2 + 9(1 (72 - 236?)E + 9(836? - fg,) 
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a{E) = 9^2 _ 8ig^^ + 27(1^2 - 29e?), 

c{E) =E^- 2161^2 + 9(1^2 - 23el)E + 9(766? - 7^3), 

Q{E) = UIoH'^'KE), ijW(E) = l, H(^\E) = E + 9e^, 

H'^^^E) = E^ + (2463 + 2762)^^ + (-14461 - 2886263 - 456^)E 
-11526i - 30246261 - 28086^63 - 71163, 

H^^\E) = E^ + (2462 + 2763)^2 + (-14461 - 2886263 - 4561)^ 

-115263 - 3024636^ - 28086^62 - 7116^, 

m(°) (E) = E^ - 5061^3 + (f ^2 + 3396^)^2 + (82086? + ^^3)^ 

9375^2 19215 „ ^ 971 c;q^4 

Ht^'^ (E) =E^- 1861^3 + ( 3|1(^2 - 21336^)^2 + (382326? + ^(73)^ 

+ 5Cp^2 + 60669^3e, _ 2178096^, 

He{E) = E-fei. 



6.5.6. The case (/q, h, I2, h) = (3, 2, 1, 0). 
S(x, E) = 225p(x)3 + 45(E - 262 - Qei)p{xf 

+3(2^2 _ 8(^3g^ ^ g^)^ _ ^2862 - 4262 + 73626i)p(x) 

+9{E + 461 - 1262)p(x + uif + 3{E + Aci - 12e2){E - lOci + 2e2)p{x + Ui) 
+ (^2 + 2(62 - 156i)E + 146261 + 5762 ^ 8962)p(a; + 102) + E^ - 7(62 + 36i)E2 

+ (-11862 + 2126162 - 156^)E + 9(26626? - 356^ + 896162 + 426?) 

a{E) = 10^2 - 4(1362 + 306i)E - 66^ - 2806? + 752626i, 

c{E) =E^- 7(62 + 36i)E2 + (23O6162 + 362 _ iooe?)E - 3816^ + 5676i62 + 3006?, 

QiE) = ULoH^'KE), H('\E) = E + 461 - 1262, H('KE) = 1, 

i7(2)(^) = ^4 _ 4(4g^ ^ g^)^3 ^ (152ei62 - 2806? - 1061)^2 

+ (16486?62 - 7406^ + 17926i6i - 3206?)E 
-32796^ + 55206? + 352626? + 2872616^ + 86966?ei, 

/J(3)(^) =E^ + 2(62 - 156i)^ + 146261 + 5762 ^ 39^2^ 
m(0)(E) = E^ - (4461 + f 62)^3 + (348g^g^ ^ 504g2 _ I^e^E^ 
^^_14016g2g^ - 74863 + 18366162 + ^6?)E 



87123^.4 87536 g4 + 64O8O626? - ^ei63 + 279606?62, 



Ht^'^\E) = E^ + (4362 - 376i)E2 + (9566? - 3306i62 - 6896^)E 
+59736^ - 17608626? + 292636i62 - 180206?, 

He{E) =E'- (f 61 + f 62)i? - ^6? + 1246261 - 2|3e2. 



6.5.7. The case (/q, h, hi h) = (3, 2, 2, 1). 
S(x, E) = 225p{xf + 45{E + Ae^)p{xf + 6{E^ + ScsE + dlej - ^g2)p{x) 

+9{E + 61 - 2462)p(a; + cui^ + 3{E + 61 - 2462)(^ - 661 + 462)p(a; + cui] 
+9(^ + 62 - 246i)p(a; + 102^ + 3{E + 62 - 246i)(E - 662 + 46i)p(a; + 102) 
+ {E + 61 - 2462)(^ - 2461 + 62)^(0; + CU3) 

+^3 + 1463^2 + (349g2 _ 5|3^^)^ _ 2|lg^^^ _ 1|7^^ 

a{E) = 13^2 ^ i4gg^^ ^ 234162 - 760(^2, QiE) = ULo H^'K^), 
c{E) =E^ + UesE^ + (3496^ - 133^2)^ - 3363^2 - 756^3, 
i7(o)(E) = l, i7«(E) = E + 6i-2462, H^^\E) = E - 24:61 + 62, 
H(^\E) = E^ + 563^^ + (-164(72 + 25061)^3 + {If 3263 - ^gs)E^ 

H-'-^92el + ^^363 + 1216(72)^ 

+ '-W9le3 + '-^93el + '-^9392, 
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+ (=^363 - '-^92el + 36480^1)^2 

,/ 5102001 ^2^ 26698161 „ „ i 337m23 „ „2\ p 

+ 1 To^fi'2^3 10^9392 H -^ — 93e3)^ 

I 3200000 ^3 51057923 ^2^2 14045111gr ^2 49594475 ^ „ ^ 
"^ 13 — 5'2 208 — ^^2^3 208 ds 10^93,92^3, 

Ht^^\E) = E^ + lOOes^^ + (-1914e^ + 37Qg2)E^ - (8519^3 + 5271^263)^ 
+ 2111393^^^2 _ 80000(^1 - ^-^^gses, 

Hd{E) = E^ + Se^E' + (195ei - 208g2)E + ^-^e.g^ + ^^3. 

6.5.8. The case {lo,li,l2,h) = (3,3,0,0). 

E{x, E) = 225(p(a;)3 + p{x + uj^f) + 45(E - I2ei){p{xf + p(a; + cji)^) 
+ (6^2 - 144eiE + ^(?2 - 4866?) (p(a:) + p(x + c^i)) 
+E^ - 4261^2 + (f ^2 - 99e2)E + 2808e? + 225^3, 
a{E) = 12(^2 - 24eiE - 81e? + 15^2), Q{E) = H'^^\E)H^^\E) 
c{E) =E^- 4261^2 + 9(5^2 - lleDE + 2448e? + 360^3, 
H^^\E) = {E- 12ei)(E2 - 36ei^ + 60^2 - 432e?), 
i^(i)(E) = E*- 36eiE^ + 66(^2 - 96^)^2 ^ (48goe3 + 324^3)E 

+225^2 ^ 6966^361 + 27297ef , 
iJ(2)(^) = 1^ ij(3)(E) = 1, 

Ht^^\E) = (^2 - 24eiE - 981e? + 75^2) (^^ - ^Ae^E + 7296? - 75^2) 
Ht^^\E) = E^ - 1861^3 + (525^2 - 6291e2)E2 + (173448ef + 2imQg^)E 
+22500^1 + 153900^361 - 2007504ef , 

He{E) = E'- ^eiE3 + (^^^ + 16M9g2)^2 + (12||88g3 + 2^10^3)^ 
22500 ^2 607500 ^ ^ 2111184 ^4 

— rr92 n— fi'3ei n — ^i- 

6.5.9. r/ie case (/q, /i, /2, ^3) = (3, 3, 1, 1). 

S(x, E) = 225{p{xf + p{x + ujif) + 45{E - 10ei)(p(x)2 + p{x + ujif) 
+ (6^2 - f ^2 - 525e2 - l20Eei){p{x) + p{x + coi)) 
+ (^2 _ 50e^^ + 25(13e2 + g2)){p{x + 102) + p{x + U3)) 
+E^ - 3561^2 + (1^2 - 200e2)E + 300(13e3 - 2^3) 
a{E) = UE^ - 340eiE - 400e2 + 80^2, 

c{E) =E^- 3561^2 + 10(^2 - 20el)E + 240(15ef - 2^3), Q{E) = H^^\E)H^^\E) 
H^'^\E) = (^2 + 4eiE - SOej - 20g2){E^ - 24eiE^ + 16(^2 - 39e2)E - 1280(ef + ^3)), 
if«(E) = E2-50eiE + 25(13e2 + ^2), H(^\E) = I, H^^\E) = l, 
Ht^^\E) = (^2 - 50eiE + 200(5e2 - g2)){E^ - fe,E + ^{g2 - 236?)) 
Ht^^\E) = E^ + 48eiE^ + 12(-796e2 + 67^2)^^ + 320(33^3 + 188el)E^ 

+9600(-^| + 27^361 + 334ef)E2 + 384000(51^36? - 64ef - 5^3^2)^ 
+640000(351^36? - 776ef + 30^| + gl), 
Hd{E) = E^- ^ei^s - i|(700e2 + U7g2)E^ + ^(-27(73 + 1640e?)E3 



f^(-5^2 _ 201^361 + 886ef)E2 + 38||0D(_37^3e2 _ I02ef + 5^3^2)^ 
+ ^52p( 159^36? - 872e6 - 66^2 ^ ^3)_ 

6.5.10. The case {lo,h,l2,h) = (3,3,2,2). 
E{x, E) = 225{p{xf + p{x + uif) + 45{E - 6ei){p{xf + p{x + uj,f) 

+ (6^2 - ^g2 - 459e2 - 72Ee,){p{x) + p{x + cu,)) 

+9{E - 36ei){p{x + 102)^ + p(x + cjg)') 

+3(^ + 6ei)(E - 36ei){p{x + 002) + p{x + u^)) 

+E^ - 2161^2 + (-72^2 - 26lel)E + 3249e? - 1917^3, 
a{E) = 18(^2 - 18eiE - 123e2 - 15^2), Q{E) = H^^\E)H^^\E) 
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c{E) =E^ - 2161^2 + (-2616? - 63^2)^ + 28536? - 1773^3 , 
i7(o)(E) = E-36ei, H^^\E) = l, H^'^\E) = l, 
H^^\E) = E^ - 6eiE^ + (-405e? - 189^2)^^ + (-4212^3 - 2A8Ael)E^ 
+ (8451^1 - 55242^361 + 69903ef)E2 
+ (-293868^36? - 275238ef + 345222^3^/2)^ 
+1426653c/3e? - 7172631e? + 1019871^| + 50625^|, 
m(o) (E) = {E^ - 3661^3 + 54(37e2 - 13^2)^^ + 108(-67e? + 107^3)^ 

+81(625^1 + 7382^361 - 13199ef))(^^ + 2861^^ - 30(^2 + 556^)^2 
+468(-61e? + 5^3)^ - 9375^1 - 179370^3ei + 476289ef), 
Ht^^\E) = E^ + ISSeiE^ + (-17253e2 + 1971c/2)^^ + (-5076^3 + 45036el)E^ 
+ (-388125^1 - 4881546^361 + 15619311et)S^ 
+ (105825204^36? - 157540950ef - 10226250^3^2)^ 
+4936453389^36? - 5402344167e? - 149900625^| + 31640625^|, 
He{E) = E^- ^e,E^ - (^ef + ^^2)^' + {^^93 + ^^el)E^ 

I ^21875^2 10529514 ^ „ , 16361919 ^4\ 7712 



,;''-gA+'-^W^e\+'^^g,g,)E 



17 H2 17 ^^Hi < 17 

134787564 „ „2 ,_ 7429158 ^5 1 12656250 , 

17 '^l "T 17 y-ib 
60220^301 ^ ^3 1 41851)5049 ^6 1287140625 ^2 , 31640625 ^3 

17 — 93^1 H jt — Ci p^ 5-3 H Yr — 92- 



6.5.11. The case {lo,h,h,h) = (3,3,3,3). 
E{x, E) = 225{p{xf + p{x + ujif + p{x + tJ2)^ + p{x + uj^f) 

+A5E{p{xy + p{x + CJi)2 + p{x + CU2)^ + Pix + CU3)2) 

+^3 - 195^2^ - 2250^3, 

a{E) = 24(^2 - 60^2), c{E) = E^ - 180^2^ - 2160^3, 

Q{E) = ijW(E) = EYll=i{E^ + '24:ekE + 720ei - 240(^2), 

iy«(E) = H^^^\E) = H^^\E) = 1, 

Ht(^\E) = {E^ - 300(^2) ULiiE' + QOckE + 300(^72 - 12ei)), 
Ht^^\E) = ^12 + 2886^^11 + 144(41t/2 - 288el)E^^ - 8640(32e^ + 17^3)^^ 
+21600(-287^| - 3744c/3efc + 9216et)^® + 12960000(-64e| + ^2^3 + IGg^eDE'^ 
+4320000(-82944e^ + 629^| - 7830^| + 63936^36^)^^ 
+777600000(280^|efe + 1024e^ + 17^1^3 - 512^36^)^^ 
+972000000(7056^1^2 - 589^| - 322560c/3eE + 62208^|e^ + 294912e|)E^ 
+233280000000(^1^3 - 1920^3e| - 1248^|e^ + 25606^ - 148^1)^^ 
+2332800000000(-18432^|et + 19^| + 2016^|efe - 36864e^° + 50688^36^ - 225g^gl)E'^ 
+279936000000000(-120^|e^ + 96^2^5 _ ^4^^ ^ 35^3^^ _ 20486^^ + 2048^3e|)E 
+466560000000000(-54^|^| + 729^| + gf) 
He{E) = E'^ - ^g2E'^ + M^(73^' + 885600glE' - 1905120000 ^^^^^7 

_43|p(i37^3 + 13986(^1)^6 + 333590400000 ^2^^^5 + i26360ooott) (2160^2^, - 71gl)E^ 
+ ™™ (1188(?g - l25glgs)E^ + ^i430720ooooooo ^2(^3 _ 27gl)E^ 

^ 25194^1000000000 ^^^^(^3 _ 27^|)^ _ 4665600oddoOOOO (^3 _ 27gl)^ 



6.5.12. The cases {lo,h,l2,h) = (2,1,1,1), (2,2,1,0), (2,2,2,0). It is shown by 
calculating functions E{x,E), a{E), c{E), Q{E), H^'\E), Ht^'\E) {i = 0,1,2,3), 
H9{E) explicitly that the monodromy for the case {lo,h,h,h) = (2,1,1,1) (resp. 
{Iq, h, I2, I3) = (2, 2, 1, 0), (2, 2, 2, 0)) is the same as the one for the case (/q, h, h, h) = 
(3, 0, 0, 0) (resp. (/q, k, k, k) = (3, 1, 0, 0), (3, 1, 1,1)). 
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For the case (/q, h, I2, h) = (2, 1,1,1) we have 
E{x, E) = 9Ep{xf + 2,E^p{x) + ELi(^' + Scfe^ + 45e| - l^g2)p{x + uok) 

+E^ - Ug^E - if (73, 
and the rest of the functions are the same as the ones in the hst for the case 
(/o,/i,/2,/3) = (3,0,0,0). 

For the case (/q, /i, h^ h) = (2, 2, 1, 0) we have 
E{x, E) = 9{E - 362 - 8es)p{xy + 3{E + 4e2 + Qe3){E - 3e2 - 863)^(3;) 

+9{E - 8e2 - 3e3)p(x + toif + 3{E + 6e2 + 4e3)(E - 8e2 - 3e3)p(a; + cui) 
+ {E^ - We^E - 32e2 + 5g2)pix + UJ2) 
+E^ - 761^2 - Ag2E - 1046?^ - ^^3 - 135ef , 
and the rest of the functions are the same as the ones in the hst for the case 
(/o,/i,/2,/3) = (3,1,0,0). 

For the case (/q, h, I2, h) = (2, 2, 2, 0) we have 

E{x, E) = 9{E - I5es)p{xf + 3{E + 6e^){E - 15e3)p{x) + E^ - '-fg2E - 540g, 

+9{E - 15e2)p(a; + cu^f + 3{E + 6e2)(E - I5e2)pix + cui) 

+9{E - I5ei)p{x + U2f + 3{E + 6ei)(E - 15ei)p(x + cus), 

and the rest of the functions are the same as the ones in the hst for the case 

(/0,/l,/2,/3) = (3,1,1,1). 

6.6. The case (7 = 4. If the genus of the related hyperelhptic curve v^ = —Q{E) is 
four and the condition Iq > h > I2 ^ h is satisfied, then we have 24 cases 

(4,0,0,0), (4,1,0,0), (4,1,1,0), (4,1,1,1), (4,2,0,0), (4,2,1,0), 

I (4,2,1,1), (4,2,2,0), (4,2,2,2), (4,3,0,0), (4,3,1,0), (4,3,2,1), 

Ifo, h, h, t3) S (4^ 3^ 3^ 2), (4, 4, 0, 0), (4, 4, 1, 1), (4, 4, 2, 2), (4, 4, 3, 3), (4, 4, 4, 4), 

(3,2,1,1), (3,2,2,0), (3,3,1,0), (3,3,2,0), (3,3,3,1), (2,2,2,1). 

If 2 > /q > /i > ^2 > ^3 then the genus is no more than four, and the genus equal to 
four case is (/q, h, I2, h) = (2, 2, 2, 1). This case is related with the case (/q, h, h, h) = 
(4,0,0,0). Now we consider the cases {lo,h,l2,h) = (4,0,0,0) and {lo,hj2,h) = 
(2,2,2,1). 

6.6.1. The case {loJukJs) = (4,0,0,0). 
E{x, E) = lW25p{xy + 1575Ep{xf + 135{E^ - fg2)p{xf 

+10(^3 _ 371^^^ _ 2A5gs)pix) + E^ - '-fg2E' - '^g,E + ^^i, 
c{E) =E'- '-fg2E' - '-^gsE + '^gl 
a{E) = 10^3 _ ^g^E - 875(?3, Q{E) = UIoH^'KE), 
H^^\E) = E^ - 52g2E - 560^3, Ht^^\E) = E^ - ^^2^ + ^^3, 
if(^)(E) = ^2 - lOckE - 35e| - 7^2, {k = 1, 2, 3) 
Ht^^\E) = E^ + 55ekE^ + |(77^2 - 5406^)^^ - 490(16ef + gs)E 

+343(-27(7i - 380g3ek + 720et), {k = 1,2,3), 
He{E) = E^- if ^2. 

For the case (/q, h, I2, h) = (1, 2, 2, 2) we have 
E{x, E)= E^- fg2E' - '-^g,E + '-^gl + {E^ - 52g2E - 560gs)p{x) 

+ ELi (9(^' - lOcfc^ - 35el - 7g2)p{x + u^f 
+3(^3 _ gg^^2 _ ^g^E _ 75g2^ ^ 28^3 - 252e|)p(a; + Uk)) , 
and the rest of the functions are the same as the ones in the list for the case 
{lo,h,h,h) = (4,0,0,0). The case {lo,h,l2,h) = (2,2,2,1) is obtained by paral- 
lel transformation x ^ x + U3 from the case (/q, h, I2, h) = (1, 2, 2, 2). For this case, 
the functions except for S(a;, E) are the same as the ones in the list. 



HEUN EQUATION IV 33 

Appendix A. Elliptic functions 

This appendix presents the definitions of and the formulas for the elhptic functions. 
The Weierstrass p-function, the Weierstrass sigma-function and the Weierstrass 
zeta-function with periods (2cji,2c(J3) are defined as follows: 

(A-1) P{z) = ^ + 2^ \{z - 2muJi - 2nuj^f ~ (2mcui + 2ncj3) V ' 



{m,n)eZxZ\{(0,0)} 

z 



a[z] = z 



n 1- 



2mui + 2nuj^ 

{m,n)eZxZ\{(0,0)} ^ 

z z^ 



■ exp ( 7^-— — ^ + 



2'mijj\ + 2na;3 2(2ma;i + 2nw^Y J 
a'(z) 



a{zy 

Setting 0J2 = —oJi — uj^ and 

(A.2) Cfc = p(c^fc), % = C(^fe) (^ = 1,2,3) 

yields the relations 
(A.3) ei + 62 + 63 = ?7i + ?72 + 773 = 0, 

^1^3 - %^i = ^3^2 - ^2^3 = ^2^1 - ^1^2 = vrv^^/2, 
P(2;) = -C'(2;), {p{z)f = 4:{p{z) - ei){p{z) - e2){p{z) - 63), 
aiz + z)aiz — z) 



Pi.z) - p{z) 



aizYaiz)'^ 



The addition formula for the Weierstrass zeta function is given as 

(A,4) C(...)-C(..-C«4f}^^ 

The periodicity of functions p{z)^ C^{z) and a{z) are as follows: 
(A.5) p{z + 2uOk) = p{z), C{z + 2uJk) = C{z) + 27]k (A; = 1,2,3), 

(j{z + t + 2ojk) a{z + t) 

a{z + 2uk) = ~or{z) exp{2r]k{z + Uk)), 7 — —z — ^— = exp{2r]kt) — — 

a[z + 2LJk) o'yz) 

The constants (72 and g^ are defined by 

(A.6) 6-2 = -4(eie2 + 6263 + 6361), 5-3 = 4616263. 

The co-sigma functions cr}^[z) {k = 1,2,3) and co-p functions pk{z) {k = 1,2,3) are 
defined by 

.. „x / s / a{z + ujk) , . (Tk{z) 

(A.7) ak{z) = exp{-r]kz) — ^, pk{z) = ^^— , 

and satisfy 

(A.8) pk{zf = p{z) - ek, {k,k' = 1,2,3) 

Pk{z + 2uk') = exp{2{r]k'UJk - rikUJk'))pk{z) = {-lY'^'"' pkiz). 
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